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Abstract 

Combined analyses of the Higgs boson production and decay rates as well as its coupling 
strengths to vector bosons and fermions are presented. The combinations include the res¬ 
ults of the analyses of the H —> yy, ZZ*, WW*, Zy, bb, tt and pp decay modes, and the 
constraints on the associated production with a pair of top quarks and on the off-shell coup¬ 
ling strengths of the Higgs boson. The results are based on the LHC proton-proton collision 
datasets, with integrated luminosities of up to 4.7 fb -1 at \[s = 7 TeV and 20.3 fb _ 1 at 
Vs = 8 TeV, recorded by the ATLAS detector in 2011 and 2012. Combining all production 
modes and decay channels, the measured signal yield, normalised to the Standard Model ex¬ 
pectation, is 1.18+^4. The observed Higgs boson production and decay rates are interpreted 
in a leading-order coupling framework, exploring a wide range of benchmark coupling mod¬ 
els both with and without assumptions on the Higgs boson width and on the Standard Model 
particle content in loop processes. The data are found to be compatible with the Standard 
Model expectations for a Higgs boson at a mass of 125.36 GeV for all models considered. 
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1. Introduction 


In 2012, the ATLAS and CMS Collaborations at the Large Hadron Collider (LHC) reported the observa¬ 
tion of a new particle at a mass of approximately 125 GeV [1,2]. The discovery made in the search for the 
Standard Model (SM) Higgs boson ( H ), is a milestone in the quest to understand electroweak symmetry 
breaking (EWSB). Within the SM, EWSB is achieved through the Brout-Englert-Higgs mechanism [3- 
8 ] which predicts the existence of a neutral scalar particle, commonly known as the Higgs boson. While 
the SM does not predict the value of its mass (mu), the production cross sections and decay branching 
ratios (BR) of the Higgs boson can be precisely calculated once the mass is known. Therefore, preci¬ 
sion measurements of the properties of the new particle are critical in ascertaining whether the newly 
discovered particle is fully responsible for EWSB and whether there are potential deviations from SM 
predictions. 

At the LHC, SM production of the Higgs boson is dominated by the gluon fusion process gg —» H (ggF), 
followed by the vector-boson fusion process qq' —» qq'H (VBF). Associated production with a W boson 
qq' —> WH ( WH ), a Z boson qqlgg —» ZH ( ZH ) or with a pair of top quarks qqlgg —> ttH ( ttH ) 
have sizeable contributions as well. The WH and ZH production processes are collectively referred to 
as the VH process. Contributions are also expected from bb —» H (bbH ) and production in association 
with a single top quark (tH). The latter proceeds through either the qb —> tHq' or gb —» WtH process. 
With the present dataset, the LHC is expected to be most sensitive to the Higgs boson decays of H —> 
yy, ZZ*, WW*, tt and bb. Together they account for approximately 88% of all decays of a SM Higgs 
boson at mu ~ 125 GeV. 

The discovery of the Higgs boson was made through analyses of the bosonic decay modes in H —* yy, 
H —» ZZ* —» 4 ■£ and H — > WW* —» (’vlv (l = e, //) events. Since the discovery, these analyses have been 
improved and updated with more data [9-11]. The H —> WW* —> tvlv analysis has been supplemented 
with a dedicated VH analysis targeting H —> WW* [12]. The ATLAS Collaboration has measured the 
Higgs boson mass from the H —> yy and H —» ZZ* —» 4 C decays to be mu = 125.36 ± 0.41 GeV [13], 
reported results in the H —> tt [14] and H —* bb [15] fermionic decay modes, and published upper limits 
on the rare decays H —> Zy [16] and H —> pp [17]. Furthermore, constraints have been set on the ttH 
production rate [18-20] and on the off-shell coupling strengths of the Higgs boson [21]. These results are 
based on the full proton-proton collision data with integrated luminosities of up to 4.7 ftr 1 at a centre-of- 
mass energy yfs -1 TeV recorded in 2011 and 20.3 fb -1 at yfs = 8 TeV recorded in 2012 by the ATLAS 
detector at the LHC. A detailed description of the ATLAS detector can be found in Ref. [22]. 

This paper presents the combined results of the analyses mentioned above. These analyses are designed 
for maximum sensitivities to SM Higgs boson production from different processes, exploiting in particular 
the differences in kinematics through categorisation of the selected events. Thus the yields of different 
Higgs boson production processes and decays can be extracted. The Higgs boson coupling strengths to 
SM vector bosons and fermions in different benchmark models are probed for the measured Higgs boson 
mass of mu = 125.36 GeV. All results are obtained assuming the Higgs boson has a small total decay 
width such that its production and decay factorise. The ATLAS Collaboration has previously published 
combined studies of Higgs boson production and decay rates [23] and of spin-parity properties [24, 25] 
using diboson final states. The results are found to be consistent with expectations from the SM Higgs 
boson. Compared with the previous publication, the current results are based on the improved analysis 
sensitivities and the addition of information from other decay modes. A similar combination has been 
published by the CMS Collaboration [26]. 
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The paper is organised as follows. Section 2 briefly summarises the individual analyses that are included 
in the combinations and Section 3 outlines the statistical method and the treatment of systematic uncer¬ 
tainties used in the combinations. In Section 4, the measured Higgs boson yields are compared with the 
SM predictions for different production processes and decay modes. In Section 5, the coupling strengths 
of the Higgs boson are tested through fits to the observed data. These studies probe possible deviations 
from the SM predictions under various assumptions, motivated in many cases by beyond-the-SM (BSM) 
physics scenarios. An upper limit on the branching ratio to invisible or undetected decay modes of the 
Higgs boson is also set. Finally, a brief summary is presented in Section 6. 


2. Input analyses to the combinations 

The inputs to the combinations are the results from the analyses of H —> yy, ZZ*, WW*, rr, bb, pp and 
Zy decay modes, and of the constraints on ttH and off-shell Higgs boson production. These analyses and 
changes made for the combinations are briefly discussed in this section. The ATLAS Collaboration has 
also performed a search for the rare H —> J/i]sy decay [27] which has the potential to constrain the Higgs 
boson coupling strength to the charm quark. However, the current result does not add sensitivity and is 
therefore omitted from the combinations. Furthermore, the inclusion of the results from direct searches 
for Higgs boson decays to invisible particles, such as those reported in Ref. [28, 29], is beyond the scope 
of the combinations presented in this paper. 

The theoretical calculations of the Higgs boson production cross sections and decay branching ratios have 
been compiled in Refs. [30-32] and are summarised in Table 1. For the ggF process, the cross section is 
computed at up to NNLO in QCD corrections [33-38] and NLO in electroweak (EW) corrections [39-41]. 
The effects of QCD soft-gluon resummations at up to NNLL [42] are also applied. These calculations are 
described in Refs. [43-47]. For the VBF process, full QCD and EW corrections up to NLO [48-50] and 
approximate NNLO [51, 52] QCD corrections are used to calculate the cross section. The cross sections 
of the WH and ZH (qq —> ZH ) are calculated including QCD corrections up to NNLO [53, 54] and EW 
corrections up to NLO [55,56] whereas the cross section of the gg —> ZH process is calculated up to NLO 
in QCD corrections [57, 58]. The cross section for ttH is computed up to NLO in QCD [59-62]. For the 
bbH process, the cross section is calculated in QCD corrections up to NLO [63-65] in the four-flavour 
scheme and up to NNLO [66] in the five-flavour scheme with the Santander matching scheme [67]. 
The cross sections of the tH processes used are calculated at up to NLO in QCD corrections [68, 69]. 
The PDF sets used in these calculations are CT10 [70, 71], MSTW2008 [72], NNPDF2.1 [73, 74] and 
NNPDF2.3 [75] following the prescription of Ref. [76]. The decay branching ratios of the Higgs boson 
are calculated using the Hdecay [77, 78] and Prophecy4f [79, 80] programs, compiled in Ref. [81]. 

All analyses use Monte Carlo (MC) samples to model the acceptances of the Higgs boson events. Table 2 
summarises the event generators and parton distribution functions (PDF) used for the analyses of the 
sfs = 8 TeV data. The modelling at \fs = 1 TeV is similar, with one notable difference of Pythia6 [83] 
replacing Pythia8 [84]. The ggF and VBF production of the Higgs boson are simulated with the next- 
to-leading order (NLO) matrix-element Powheg program [85-89] interfaced to either Pythia6 or Pythia8 
for the simulation of the underlying event, parton showering and hadronisation (referred to as the shower¬ 
ing program). The Higgs boson transverse momentum distribution from ggF production is reweighted 
to match the calculation of HRes2.1 [90, 91], which includes QCD corrections up to the next-to-next- 
to-leading order (NNLO) and next-to-next-to-leading logarithm (NNLL) in perturbative expansions. Fur¬ 
thermore, ggF events with two or more jets are reweighted to match the transverse momentum distribution 
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Table 1: SM predictions of the Higgs boson production cross sections and decay branching ratios and their uncer¬ 
tainties for m H = 125.36 GeV, obtained by linear interpolations from those at 125.3 and 125.4 GeV from Ref. [32] 
except for the tH production cross section which is obtained from Refs. [20, 82]. The uncertainties of the cross 
sections are the sum in quadrature of the uncertainties resulting from variations of QCD scales, parton distribution 
functions and a s . The uncertainty on the tH cross section is calculated following the procedure in Refs. [20, 32], 


Production 

Cross section [pb] 

Decay channel 

Branching ratio [%] 

process 

yfs = 7 TeV 

yfs = 8 TeV 

H —> bb 

57.1 ± 1.9 

ggF 

15.0+ 1.6 

19.2 + 2.0 

H -» WW* 

22.0 ± 0.9 

VBF 

1.22 + 0.03 

1.57 + 0.04 

H -» gg 

8.53 + 0.85 

WH 

0.573 + 0.016 

0.698 + 0.018 

H —» TT 

6.26 ± 0.35 

ZH 

0.332 + 0.013 

0.412 + 0.013 

H -» cc 

2.88 + 0.35 

bbH 

0.155 + 0.021 

0.202 ± 0.028 

H -» ZZ* 

2.73 + 0.11 

ttH 

0.086 ± 0.009 

0.128 + 0.014 

H —* yy 

0.228 + 0.011 

tH 

0.012 + 0.001 

0.018 + 0.001 

H ->Zy 

0.157 + 0.014 

Total 

17.4+ 1.6 

22.3 ± 2.0 

H —»ju/r 

0.022 ± 0.001 


from MiNLO HJJ predictions [92]. The WH and ZH (qq —> ZH) production processes are simulated with 
the leading-order (LO) Pythia8 program. The gg —> ZH process contributes approximately 8% to the 
total ZH production cross section in the SM. For most of the analyses, the process is modelled using 
qq —> ZH of Pythia8. Only the VH analysis in the H —* bb decay mode specifically models gg —> ZH 
production using Powheg [85-87] interfaced to Pythia8. The ttH process is modelled using the NLO 
calculation in the FIELAC-Oneloop package [93] interfaced to Powheg and Pythia8 for the subsequent 
simulation. The tH production process is simulated using MadGraph [94] interfaced to Pythia8 for 
qb —> tHq' and using MadGraph5_aMC@NLO [82] interfaced to Herwig++ [95] for gb —> WtH. The 
bbH production process contributes approximately 1 % [96] to the total Higgs boson cross section in the 
SM. It is simulated with the MadGraph5_aMC@NLO program for some analyses. The event kinematics 
of ggF and bbH production are found to be similar for analysis categories that are most important for 
bbH. Thus the acceptance times efficiency for bbH is assumed to be the same as for ggF for all analyses. 
The PDF sets used in the event generations are CT10 [70] and CTEQ6L1 [97]. All Higgs boson decays 
are simulated by the showering programs. 

Throughout this paper, the signal-strength parameter p is defined as the ratio of the measured Higgs boson 
yield to its SM expectation: 


^ = 


cr x BR 
(cr x BR)sm 


( 1 ) 


Here cr is the production cross section of the Higgs boson. For a specific production process i and decay 
channel /, i.e., i —> H —> /, the signal-strength parameter is labelled as and can be factorised in terms 
of the signal strengths of production (//,■) and decay (fif). 


= 


0~i X B R y 


= HiX/Uf, with m = 


<Ti 


and /if ~ 


BR 


/ 


( 2 ) 


(cr i X BR/)sm (O' (Ism (BR/)sm 

Thus for each analysis category (c) as discussed later in this section, the number of signal events ( n c s ) can 
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Table 2: Summary of event generators, showering programs and PDF sets used to model the Higgs boson production 
and decays at = 8 TeV. 


Production 

process 

Event 

generator 

Showering 

program 

PDF 

set 

ggF 

Powheg 

Pythia6/Pythia8 

CT10 

VBF 

PoWHEG 

Pythia6/Pythia8 

CT10 

WH 

Pythia8 

Pythia8 

CTEQ6L1 

ZH : qcj -> ZH 

Pythia8 

Pythia8 

CTEQ6L1 

ZH : gq -> ZH 

Powheg 

Pythia8 

CT10 

ttH 

Powheg 

Pythia8 

CT10 

bbH 

MadGraph5_aMC @NLO 

Herwig++ 

CT10 

tH : ql) —> tHq' 

MadGraph 

Pythia8 

CT10 

tH : gb -> WtH 

MadGraph5 aMC @NLO 

Herwig++ 

CT10 


be written as: 


n c s = 2 2 x iu :/(BR/) S m x A c if x e c if x £ c (3) 

«' / 

where the indices i and / indicate the production processes and decays contributing to the category, Ay 
represents the detector acceptance derived from simulation of the SM process, sy is the reconstruction 
efficiency within the acceptance and the integrated luminosity for the given category c of the given 
channel. 

However, the experimental data do not allow to separately determine and /if for any given process 
since only their - product is measured. All combined fits of signal strengths presented in this paper make 
assumptions about the relationship between p, of different production processes or similarly between p/ 
of different decay modes. Thus the meaning of the signal strength depends on the assumptions made. Nev¬ 
ertheless, the production and decays can be factorised using the ratios of cross sections and of branching 
ratios as discussed in Section 4.4. 

Leptons (t) refer to electrons or muons unless specified otherwise; the symbols ri ep and Th a d refer to 
t leptons identified through their decays to leptons or hadrons; and variables p T , E \ and f?™ ss refer 
to transverse momentum, transverse energy and missing transverse momentum, respectively. Notations 
indicating particle charges or antiparticles are generally omitted. 

The ATLAS experiment uses a right-handed coordinate system with its origin at the nominal interaction 
point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the 
IP to the centre of the LHC ring, and the y -axis points upward. Cylindrical coordinates (r, (p) are used in 
the transverse plane, (b being the azimuthal angle around the beam pipe. The pseudorapidity is defined in 
terms of the polar angle 6 as ij = - In tan(0/2). 

Table 3 gives an overview of the analyses that are inputs to the combinations and their main results, as 
published. An essential feature of these analyses is the extensive application of exclusive categorisation, 
i.e., classifying candidate events based on the expected kinematics of the different Higgs boson production 
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processes. The categorisation not only improves the analysis sensitivity, but also allows for the discrim¬ 
ination among different production processes. Figure 1 summarises the signal-strength measurements of 
different production processes that are used as inputs to the combinations. 


ATLAS 

Individual analysis 


H —» yy 


Overall: p = 1. 
ggF: p = 1. 
VBF: p = 
WH: p = 
ZH: p = 


1?+ 0.27 

-0.27 

32 +0 ' 38 

-0.38 

0 . 8 m7 

-0.7 

i.o* 16 
-1.6 
n h +3.7 


m„ (GeV) 


Input measurements 

± la on p 


125.4 

125.4 

125.4 

125.4 

125.4 


H -»ZZ* 


Overall: p = 1.44 


ggF+ttH: p = 
VBF+VH: p = 


1 .7™' 

-0.4 

0.3 +ie 


125.36 

125.36 

125.36 


H -> WW* 


Overall: p = 1.16 


ggF: p = 0. 
VBF: p = 1 
VH: p = 


i-98 


- 0.24 

-0.21 

-0.29 

0.26 

2g+0.55 

-0.47 

3.0 +1 ' 6 

-1.3 


125.36 

125.36 

125.36 

125.36 


Overall: p = 1 
ggF:p = 
VBF+VH: p = 1. 


.43 


-0.43 

-0.37 

2.0 +1 ' 5 

-1.2 

2 4 +o-59 

-0.54 


125.36 

125.36 

125.36 


VH -> Vbb 


Overall: p = 0.52 +a40 

-0.40 
^ ^+0.65 
-0.61 
! 05*0-52 
-0.49 


WH: p = 1." 
ZH: p = O.C 


125.36 

125 

125 




Overall: p = -0.7 +d 


125.5 


H ->Zy 


Overall: p = 2.7 


125.5 


ttH 


bE>: p = 1.5*^ 
Multilepton: p = 2.1 + ’ 4 
Y7:p = 1.3 + f“ 


125 

125 

125.4 


i. . I 


is = 7 TeV, 4.5-4.7 fb' 1 
is = 8 TeV, 20.3 fb' 1 


-2 0 2 4 

Signal strength (p) 


Figure 1: Summary of the signal-strength measurements, as published, from individual analyses that are inputs 
to the combinations. The Higgs boson mass column indicates the inn value at which the result is quoted. The 
overall signal strength of each analysis (black) is the combined result of the measurements for different production 
processes (blue) assuming SM values for their cross-section ratios. The error bars represent ± lcr total uncertainties, 
combining statistical and systematic contributions. The green shaded bands indicate the uncertainty on the overall 
signal strength obtained by each analysis. The combined signal strength of the H —> yy analysis also includes the 
ttH contribution which is listed separately under ttH production. 
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Table 3: Overview of the individual analyses that are included in the combinations described in this paper. The 
signal strengths, the statistical significances of a Higgs boson signal, or the 95% CL upper limits on the Higgs 
boson production rates or properties are also shown wherever appropriate. A range is quoted for the upper limit 
on the off-shell signal strength, depending on the assumption for the continuum gg —» WW/ZZ cross section. 
These results are taken directly from the individual publications. Results of the on-shell analyses are quoted for 
nin = 125.36 GeV except that m H = 125.5 GeV is assumed for the H — * Zy and // —> pp analyses and that 
nin = 125 GeV is used for the ttH searches with H —> bb and ttH —> multileptons. The luminosity used for the 
= 7 TeV VII (—> bb) analysis differs slightly from the values used for other analyses because a previous version 
of the luminosity calibration was applied. The significance is given in units of standard deviations (s.d.). The 
numbers in parentheses are the expected values for the SM Higgs boson. The ttH analysis in the H —> yy decay is 
part of the H —> yy analysis. It is included separately under the ttH production for completeness. The checkmark 
(/) indicates whether the analysis is performed for the respective yfs - 7 and 8 TeV dataset. 


Analysis 

Signal 

f £dt [ftr 1 ] 

Categorisation or final states 

Strength p 

Significance [s.d.] 

7 TeV 

8 TeV 

H — > yy [9] 

1.17 + 0.27 

5.2 (4.6) 

4.5 

20.3 

ttH: leptonic, hadronic 



/ 

/ 

VH: one-lepton, dilepton, LT" SS , hadronic 


/ 

/ 

VBF: tight, loose 



/ 

/ 

ggF: 4 p-\\ categories 



/ 

/ 

H — > ZZ* — > M [10] 

i 44+0.40 
i, ^ + -0.33 

8.1 (6.2) 

4.5 

20.3 

VBF 


/ 

/ 

VH: hadronic, leptonic 



/ 

/ 

ggF 



/ 

/ 

//-> Wr [11, 12] 

1 16 +0 - 24 
a ' au -0.21 

6.5 (5.9) 

4.5 

20.3 

ggF: (0-jet, 1-jet) ® (ee + pp, ep) 


/ 

/ 

ggF: > 2-jet and ep 




/ 

VBF: > 2-jet ® (ee + pp, ep) 



/ 

/ 

VH: opposite-charge dilepton, three-lepton, four-lepton 


/ 

/ 

VH: same-charge dilepton 




/ 

H -> TT [14] 

1 4^+0.43 

1 '^ J -Q37 

4.5 (3.4) 

4.5 

20.3 

Boosted. Ti e pTi e p, Ti e pThacb ^"had^"had 



/ 

/ 

VBF. Ti e pTi e p, TjepT^ad? ^"had^"had 



/ 

/ 

VH -» Vbb [15] 

0.52 ± 0.40 

1.4 (2.6) 

4.7 

20.3 

OH (ZH -> vvbb): /Vj et = 2,3, A btag = 

= 1,2, p\ e 100-120 and > 120 GeV 

/ 

/ 

U ( WH (vbb): V jet = 2,3, A bt ag = 

= 1,2 ,p\ < and >120 GeV 

/ 

/ 

2( (ZH abb): TVjet = 2,3, V btag = 

: 1,2, pV < and >120 GeV 

/ 

/ 

95% CL limit 

H-^Zy [16] 


P< 1K9) 

4.5 

20.3 

10 categories based on A?/ Zy and p Tt 



/ 

/ 

H^pp [17] 


p < 7.0 (7.2) 

4.5 

20.3 

VBF and 6 other categories based on J] fl and p'^ 


/ 

/ 

ttH production [18-20] 



4.5 

20.3 

H —> bb: single-lepton, dilepton 


p < 3.4 (2.2) 


/ 

ttH —>multileptons: categories on lepton multiplicity 

p < 4.7 (2.4) 


/ 

H —> yy : leptonic, hadronic 


p < 6.7 (4.9) 

/ 

/ 

Off-shell H* production [21] 

I 1 

<5.1-8.6 (6.7-11.0) 


20.3 

//* —> ZZ —> 4( 




/ 

//* —> ZZ —> 2^2v 




/ 

H* WW evpv 




/ 
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2.1. H —» 77 


In the H —» 77 analysis, described in detail in Ref. [9], the Higgs boson signal is measured in events with 
at least two isolated and well-identified photon candidates. The leading and subleading photon candidates 
are required to have Ej/m yy > 0.35 and 0.25, respectively, where m yy is the invariant mass of the two 
selected photons. The diphoton candidate events are grouped into twelve exclusive categories separately 
for the sfs = 1 and 8 TeV datasets; the order of categorisation is chosen to give precedence to production 
modes with the most distinct signatures. Each category is optimised by adjusting the event selection 
criteria to minimise the expected uncertainty on the signal yield of the targeted production mode. 

The first two categories are designed for ttH production based on the topology of leptonic and hadronic 
decays of the associated ft pair. They are described in Section 2.8. The next four categories are optimised 
for VH production, targeting one-lepton, dilepton, f?™ ss , and hadronic signatures of W and Z boson 
decays. Events from VBF production arc identified by requiring two well-separated and high-px jets and 
little hadronic activity between them. A boosted decision tree (BDT) [98, 99] algorithm is employed to 
maximise the VBF signal and background separation. Events are sorted into two categories with different 
VBF purities according to the output value of the BDT. Finally, the remaining events are separated into 
four categories based on the pseudorapidities of the photons and the pxt of the diphoton system [9], the 
diphoton momentum transverse to its thrust axis in the transverse plane. 

For most of the categories, the background is composed of a mixture of y-jet and jet-jet events, where one 
or two jets are misidentified as photons, and yy events. In particular the yy background is dominant and 
irreducible. The Higgs boson signal is extracted from maximum-likelihood fits of a narrow resonance 
plus continuum background models to unbinned diphoton invariant-mass distributions observed in the 
different event categories. In the fit, the signal is modelled by the sum of a Crystal Ball function [100] 
and a smaller but wider Gaussian component while the backgrounds are modelled by category-dependent 
exponential functions of first- or second-order polynomials. 


2.2. H -> ZZ* -> At 

The H —> ZZ* —> At analysis, described in detail in Ref. [10], has a high signal-to-background ratio, 
which is about two for each of the four final states considered: 4p, 2e2fi, 2/j2t\ and Ae, where the first 
lepton pair has an invariant mass closer to the Z boson mass. The analysis selects Higgs boson candidates 
by requiring two pairs of isolated, same-flavour and opposite-charge leptons with one of the two pairs 
having a dilepton invariant mass in the range 50-106 GeV. 

To measure the rates of different production processes, each H —> ZZ* —> At candidate is assigned to one 
of four categories depending on event characteristics beyond the four selected leptons. The VBF category 
consists of candidates with two additional jets with dijet mass nijj >130 GeV. The events failing this 
selection are considered for the VH -hadronic category, where the dijet mass is required to be 40 GeV < 
rrijj < 130 GeV. Events failing the VBF and V 7 //-hadronic categorisation criteria are considered for 
the V 7 //-leptonic category with the requirement of an additional lepton. Finally, the remaining events 
are assigned to the ggF category. The separation of VBF and VH production from the dominant ggF 
production mode is improved by exploiting two BDT discriminants trained on the jet kinematics, one 
for the VBF category and the other for the VH -hadronic category. A third BDT discriminant based 
on the four-lepton kinematics is used to improve the separation between the ggF signal and its main 
background. 
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The largest background comes from continuum ZZ* production and is estimated using simulation norm¬ 
alised to the SM next-to-leading-order cross-section calculation. For the four-lepton events with an in¬ 
variant mass, rriM, below about 160 GeV, there are also important background contributions from Z+jets 
and tt production with two prompt leptons, where the additional charged lepton candidates arise from 
decays of hadrons with b- or c-quark content, from photon conversions or from misidentified jets. Their 
contributions are estimated with data-driven methods. 

For each category, the signal is extracted from a maximum-likelihood fit to either the distribution 
(VH categories) or the combined two-dimensional distributions of and a BDT discriminant (ggF and 
VBF categories). The four-lepton mass range of 110 GeV < < 140 GeV is included in the fits. 


2.3. H -> WW* 

Analyses targeting the ggF, VBF, and VH production modes [11, 12] are performed for the H —> WW* 
decay channel. The ggF and VBF production processes are explored through the H —> WW* —> tvlv 
decay and the VH process is studied in final states with two or more leptons. 

The analysis of the ggF and VBF production processes [11] selects the signal candidate events by re¬ 
quiring two oppositely charged leptons. Candidates are categorised according to the number of jets (Vj cl ) 
and to the flavours of the leptons. The 7Vj e t categorisation separates the large top-quark production back¬ 
ground from the ggF signal while the categorisation by lepton flavours isolates the challenging Drell-Yan 
background in the same-flavour categories. The categories targeting ggF production include Vj e t = 0, 1 
and > 2 and are further divided into the same- and different-flavour leptons for Vj et = 0, 1. Only the 
different-flavour leptons are considered for Vj e t > 2. The categories targeting VBF production require 
Ajet > 2, separately for the same- or different-flavour leptons. The primary background processes are 
WW, top quark (tt and Wt), VT+jets, Drell-Yan, and other diboson (WZ, Wy, Wy* , and ZZ) production. 
Most of the background contributions are estimated using data. For the ggF categories, the final signal 
region is selected by requiring the dilepton mass mu <55 GeV and their azimuthal angular separation 
A <pet <1.8 and the signal is extracted through a combined fit to the transverse mass distributions of the 
dilepton plus E' | n,ss system in both the signal and control regions of different categories and lepton fla¬ 
vours. For the VBF categories, a BDT combining information such as rapidity separation and mass of 
the two leading jets and the dilepton angular separation, is used as the final discriminant, from which the 
signal is extracted. 

The VH analysis [12] is optimised for different lepton multiplicities: opposite-charge dileptons, same- 
charge dileptons, three and four leptons. Most final states are required to have Z?™ ss and events with a 
^-tagged jet are vetoed. Dilepton final states target VH production with the H —> WW* decay with two 
bosons decaying leptonically and the other hadronically. The opposite-charge dilepton final state selects 
events with two or more jets, with the value of m-y required to be close to the W and Z boson masses. 
Similar to the ggF Vj e t > 2 category, the dominant background is from top quark production. The same- 
charge dilepton category accepts events with either one or two jets. The dominant backgrounds are from 
WZ, Wy ( *\ and IT+jets production. The three-lepton final state targets WH with H —» WW* and has the 
highest sensitivity of the four final states. The three leptons are required to have a net charge of ±1 and 
the event can have at most one jet. The dominant background process is WZ production and is reduced 
with a Z —» It veto. The four-lepton category is designed to accept events from ZH production with the 
H —> WW* decay. The net charge of the leptons is required to be zero and at least one pair of leptons is 
required to have the same flavour, opposite charge, and an invariant mass close to the Z boson mass. The 
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dominant background is SM ZZ* production. In the three-lepton category, the signal yield is extracted 
through fits to distributions of a BDT or the minimum separation in the ij - <p plane between opposite- 
charge leptons depending on the lepton flavours. For other categories, the event yields are used, without 
exploiting information on the shapes of distributions. 


2.4. H — » tt 

The H —» tt analysis [14] considers both the leptonic (ri ep ) and hadronic (jhad) decays of the r lepton. 
Three sub-channels (Ti ep Ti ep , ri ep Th a d and Th a dTh a d) are defined by orthogonal requirements on the number 
of reconstructed hadronic r decays and leptons (electrons or muons) in the event. 1 

Candidate events are divided into boosted and VBF categories. The boosted category targets signal events 
where the Higgs boson is produced with a large boost, primarily from the gluon fusion process, and 
requires the transverse momentum of the reconstructed Higgs boson candidate to be greater than 100 GeV. 
The VBF category contains events with two jets separated in pseudorapidity and targets signal events 
produced through the vector boson fusion process. A separate BDT is then employed in each category 
and sub-channel to discriminate signal from background, utilising between five and nine input variables, 
chosen in order to exploit discriminating features such as Higgs boson decay properties, event activity, and 
the VBF topology in the corresponding category. One of the most important input variables is the mass 
of the rr system, which is quite challenging to reconstruct due to the presence of at least two neutrinos in 
the final state; the Missing Mass Calculator [101] is used for this puipose. 

In all three sub-channels, the most important backgrounds are irreducible Z — > tt events, and events with 
one or two jets misidentified as t lepton decay products (primarily from multijet and Il'+jcts production). 
To estimate the Z —> tt background the embedding technique [102] is used, where Z —» qp events are 
selected in data and the reconstructed muons are replaced by simulated t lepton decays. Fully data-driven 
techniques are used for the estimation of backgrounds from misidentified t decay products, while Monte 
Carlo simulation corrected to data is used for other backgrounds, such as the top quark and Z —» It 
production. 

The signal is extracted by fitting the shape of the BDT discriminant with signal and background templates 
simultaneously in all signal regions. The fit also includes dedicated control regions enriched with top 
quark, Z —> tt and multijet events. These control regions are used to constrain normalisations of the 
corresponding backgrounds. 


2.5. VH with H -> bb 

The H —» bb decay mode is predicted in the SM to have the largest branching ratio (see Table 1). In spite 
of this large branching ratio, an inclusive search for H —» bb is not feasible because of the overwhelming 
background from multijet production. Associated production of a Higgs boson with a vector boson V 
(IT or Z), offers a viable alternative because leptonic decays of the vector boson, W —* tv, Z —> tt, and 
Z —> vv, can be efficiently used for triggering and background reduction. 

The search for associated VH production with H —> bb [15] is performed for events containing zero, 
one, or two charged leptons. Contributions from W —> tv and Z —» tt decays in which the t leptons 

1 For events with two leptons, a requirement on the invariant mass of the rr system reconstructed via the collinear approximation 
also ensures orthogonality with the H —» WW* —> tv tv analysis. 
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subsequently decay to electrons or muons are also included. A ^-tagging algorithm is used to identify 
jets from H —» bb decays. To improve the sensitivity, the three channels are each split into categories 
according to the vector-boson transverse momentum, p^, the number of jets, and the number and quality 
of the ^-tagged jets. Topological and kinematic selection criteria are applied within each of the resulting 
categories. The categories providing most of the sensitivity are those requiring two /Magged jets and 
large p x r The categories with low sensitivity are used to constrain the contributions of the dominant 
background processes. 

A binned profile maximum-likelihood fit to all categories simultaneously is used to extract the signal 
yield and the background normalisations. The most significant background sources are V +hcavy-!lavour- 
jet production and tt production. The normalisations of these backgrounds are fully determined by the 
likelihood fit. Other significant background sources are single-top-quark and diboson ( WZ and ZZ) pro¬ 
duction, with normalisations from theory, as well as multijet events. The shapes of ah backgrounds are 
estimated from simulation, except for the multijet background for which the shape and normalisation are 
obtained using multijet-enriched control samples. 

Two versions of the analysis are performed. In the dijet-mass analysis, the mass of the dijet system of b- 
tagged jets is the final discriminating variable used in the statistical analysis. In the multivariate analysis 
(MVA), which incorporates various kinematic variables in addition to the dijet mass and the /Magging 
information, the outputs of boosted decision trees provide the final discriminating variable. Since the 
MVA has higher expected sensitivity, it is chosen as the nominal analysis for the = 8 TeV dataset to 
extract the final results. For the yfs = 7 TeV dataset, only a dijet-mass analysis is performed. 

The yfs = 1 TeV VH(—> bb) analysis uses a previous version of the luminosity calibration and therefore 
has a slightly different luminosity value compared with those quoted for other analyses. However, this 
small difference is expected to have negligible effects on the results presented in this paper. 


2.6. H -> Zy 

The H —> Zy analysis [16] with Z —> It searches for a narrow peak in the reconstructed tty invariant-mass 
distribution around 125 GeV over a smooth background. The Z+y production, Z —> tty radiative decays 
and Z+jets events where a jet is misidentified as a photon dominate the background contributions. 

The analysis selects two isolated leptons of same flavour and opposite charge and one isolated photon. 
Due to the kinematics of the decay, low pj thresholds are applied to the leptons and the photon. The 
invariant mass of the dilepton system must satisfy met > m z - 10 GeV and the three-body invariant mass 
must be consistent with the mass of the Higgs boson. To enhance the sensitivity of the analysis, events are 
classified into categories with different signal-to-background ratios and invariant-mass resolutions, based 
on the pseudorapidity difference A pzy between the photon and the Z boson and pj t , the component of the 
Higgs boson candidate pj that is orthogonal to the Zy thrust axis in the transverse plane. 

The final discrimination between signal and background events is based on a simultaneous likelihood fit 
to the mtty spectra in each category, separately for the \fs = 7 and 8 TeV datasets. Similar to the H —> yy 
analysis (Section 2.1), the signal is modelled with the sum of a Crystal Ball function and a smaller 
but wider Gaussian component while the backgrounds are modelled with polynomials, or exponentiated 
polynomials depending on categories. 
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2.7. H -* nn 


The H —» pp analysis [17] searches for a naiTow peak in the dimuon invariant mass m m distribution over 
a smooth background, where the width of the signal is dominated by the experimental resolution. The 
mass spectrum is dominated by the continuously falling background due to Z/y* production, with smaller 
contributions from top quark and diboson production. 

The selected events containing a pair of oppositely charged muons are separated into seven mutually 
exclusive categories based on the VBF dijet signature, the muon pseudorapidity q p , and the transverse 
momentum of the dimuon system . The events with two or more jets that match selections designed 
for the VBF process are accepted in the VBF signal region. All other selected events are split up into six 
categories based on the values of q p and p'j. This categorisation takes advantage of the higher momentum 
resolution for muons reconstructed in the central part of the detector, and high for the expected SM 
signal. 

The m pM distribution in the 110-160 GeV region is fitted with an analytic signal-plus-background model 
separately for the yfs = 7 and 8 TeV datasets, setting a limit on the dimuon decay of the SM ffiggs 
boson with a mass of 125.5 GeV. In the fit, the signal is modelled as the sum of a Crystal Ball func¬ 
tion and a Gaussian function in all regions while the backgrounds are modelled using exponentials or 
polynomials. 


2.8. ^//production 

Searches for qqlgg —> ttH production have been performed with three analyses targeting the Higgs boson 
decays H —» bb, H —» {WW*, rr, ZZ*) —» leptons, and H —» yy. The search in the H —> yy decay mode 
uses both yfs = 1 and 8 TeV data, while the other two use only the = 8 TeV data. 

The search for ttH production with H —» bb [18] considers two separate selections optimised for single¬ 
lepton and dilepton final states of tt decays. In the single-lepton channel, events are required to have 
one isolated electron or muon and at least four jets. In the dilepton channel, events are required to have 
two opposite-charge leptons ( ee, pp or ep) and at least two jets; events consistent with originating from 
a Z —* tt decay are rejected. In both cases at least two ^-tagged jets are required. Candidate events are 
categorised according to the jet and b- jet multiplicities with a total of nine (six) categories for the single¬ 
lepton (dilepton) final states. The background is dominated by tf+jets events, with increasing fractions of 
ttbb and ttcc at the higher b- jet multiplicities characteristic of signal events. The analysis uses a neural 
network to discriminate signal from background in the most signal-like categories. Simpler kinematic 
discriminants are used in background-like categories. 

The ttH search with H —* WW* , tt and ZZ* decays [19] exploits several multilepton signatures resulting 
from leptonic decays of vector bosons and/or the presence of r leptons. The events are categorised by 
the number of reconstructed electrons or muons and hadronic r candidates. The five channels used in 
this combination are: one lepton with two hadronic r candidates, two same-charge leptons with zero or 
one hadronic r candidate, three leptons, and four leptons. The largest backgrounds in the analysis are 
non-prompt leptons, primarily arising from semileptonic 6 -hadron decays in tt events; electron charge 
misreconstruction in events where opposite-charge leptons are produced; and the production of ttW and 
ttZ (CtV). The potential signal is determined from the numbers of observed events in data and of the 
estimated number of background events. 
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The ttH search in the H —> yy channel [20] is part of the H —* yy analysis (see Section 2.1) and employs 
the same diphoton selection. The leptonic as well as fully hadronic decay signatures of the tt system are 
considered. The leptonic selection requires at least one lepton and one ^-tagged jet as well as .E™ ss . In 
the hadronic selection, different combinations of jet and ^-tagging multiplicities are applied to improve 
the signal sensitivity. The small contribution from ggF, VBF and VH production is estimated from Monte 
Carlo simulation. The ttH signal is extracted from a fit to the observed diphoton mass distribution. 


2.9. Off-shell Higgs boson production 

Measurements of the H* —> ZZ and H* —> WW final states in the mass range above the 2my and 2 rriw 
thresholds (off-shell region) provide a unique opportunity to measure the off-shell coupling strengths of 
the observed Fliggs boson, as discussed in Refs. [103-106]. The ZZ — * At, ZZ —> 2C2v and WW —> ev/j.v 
final states in the = 8 TeV dataset are used in these measurements, detailed in Ref. [21]. Assuming the 
relevant Higgs boson coupling strengths arc independent of the energy scale of Higgs boson production, 
a combination with the on-shell measurements can be interpreted as a constraint on the total width of the 
Higgs boson. 

The analysis in the ZZ —> At final state follows closely the Higgs boson measurements in the same 
final state, described in Section 2.2, with the same object definitions, event selections and background 
estimation methods. The off-peak region is defined to include the range 220 GeV < m^e < 1000 GeV. 
Like the H — » ZZ* —» At analysis, the background is dominated by qqlgg —» ZZ production. A matrix- 
element-based discriminant [21] is constructed to enhance the gg —* H* —> ZZ signal and is used in a 
binned maximum-likelihood fit for the final result. 

The analysis in the ZZ —> 2(2v channel follows closely the ZH analysis with the Higgs boson decay¬ 
ing to weakly interacting particles [28], with the same object definitions. As the analysis is performed 
inclusively in the number of jets in the final states, kinematic cuts are optimised accordingly. SM ZZ 
and WZ production are the major backgrounds. The transverse mass (n?| z ) [21], reconstructed from the 
momentum of the dilepton system and the missing transverse momentum, is chosen as the discriminating 
variable. Events in the range of 380 GeV < nrp < 1000 GeV are used in a binned maximum likelihood 
fit for the final result. 

The analysis in the WW —> ev/uv channel follows closely the Higgs boson measurements in the oppositely 
charged electron-muon pair final state, described in Section 2.3, with the same object definitions. The 
analysis is performed inclusively in the number of jets in the final state, and selections are optimised for 
the off-shell region with revised background estimation methods. Top quark pairs and WW events con¬ 
stitute the major backgrounds. In order to isolate the off-shell Higgs boson production while minimising 
sensitivity to higher-order QCD effects on gg —> WW kinematics, a new variable R% [12], defined as the 
weighted combination of the dilepton mass and the transverse mass of the dilepton and £™ ss system, is 
constructed to select the signal region. The final results are obtained from the numbers of events observed 
in the data and expected from background processes in the signal region of Rg > 450 GeV. 


2.10. Modifications of analyses 

To ensure a consistent interpretation of all inputs in terms of Higgs boson coupling strengths, several 
minor modifications were made to the inputs of these combinations with respect to their previously pub¬ 
lished versions: 
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• The upper limits on the H —> Zy and H —» pp decays and the results of the ttH searches in H —> bb 
and ttH — * multilepton decays have been updated to assume a Higgs boson mass of 125.36 GeV. 

• In some individual analyses, cross-feed of other Higgs boson decays occurs: in the VH, H —> WW* 
selection cross-feed of H —» rr and H —> ZZ* occurs (whereas this cross-feed is negligible in the 
ggF and VBF H —> WW* analyses where a veto on the reconstructed rr mass is applied). Similarly, 
there is cross-feed from H —> WW* in the H —» rr analysis. In such cases, this cross-feed was 
treated as background in the relevant individual channel analyses. For the combinations described 
in this paper, such events are interpreted as signal from the corresponding Higgs boson decay. 

• The rate of gg —> ZH events in the VH channels is parameterised in terms of Higgs boson coupling 
strengths to Z bosons and top quarks, following the calculations of Ref. [58] for yfs = 7 and 8 TeV. 

• The rate of tH events in all the ttH channels is parameterised in terms of Higgs boson coupling 
strengths to W bosons and top quarks. 

• In the standalone analysis of the ttH channels, small contributions from Higgs boson decays to 
the cc and gg final states are explicitly modelled. To avoid spurious sensitivity due to these very 
small components in the combined analyses presented in this paper, both aforementioned decays 
are treated like H —> bb in the fits for the Higgs boson signal strength. In fits for Higgs boson 
coupling strengths, it is assumed that the coupling strengths of the H —» cc and H —> gg decays 
scale as the tt —» H and gg —* H couplings, respectively. 

• Theoretical uncertainties from QCD scales in Higgs boson signal processes have been updated to be 
consistent with the latest recommendations [32] for H —» WW*, bb, tt and Zy. No modifications 
were needed for the H—> yy and H —» ZZ* channels. 

• In channels where bbH production was not explicitly modelled, the signal strength of ggF is re¬ 
defined to include this process. In channels where bbH was modelled explicitly (H —> yy,ZZ*), 
ggF and bbH production are correlated with their ratio fixed to the SM value, allowing a consistent 
treatment of bbH production across all channels. The impact of this average scaling on the results 
is negligible since, as can be seen in Table 3, the bbH production process has a cross section which 
is only 1 % of the ggF production in the SM. 

• The off-shell analysis depends on the unknown A"-factor from higher-order QCD corrections for 
the gg —> VV background process. In the case of the very similar Higgs boson signal gg —> H* —» 
VV production process, a A'-factor between 0.5 and 2 is expected, as discussed in Ref. [21]. The 
results are given as a function of the unknown ratio of the A'-factors for gg —> VV background and 
gg —> H* —> VV signal, R^,. The range 0.5-2.0 is chosen as a systematic uncertainty on R 1 ^,. 

3. Statistical procedure 

The statistical treatment of the data is described in Refs. [107-111]. Hypothesis testing and confidence 
intervals arc based on the A(o') profile likelihood ratio [112] test statistic. The test statistic depends on 
one or more parameters of interest a, such as the Higgs boson signal strength p normalised to the SM 
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expectation (Eq. (1)), Higgs boson mass m#, coupling strength scale factors k and their ratios A, as well 
as on additional parameters 6 that arc not of interest, 


L(a , 6(a)) 
L(a, 0) 


(4) 


The likelihood functions in the numerator and denominator of the above equation are built using sums of 
signal and background probability density functions (pdfs) of the discriminating variables, introduced in 
Section 2. The pdfs are derived from MC simulation for the signal and from both data and simulation for 
the background. Likelihood fits to the observed data are done for the parameters of interest. The single 
circumflex in Eq. (4) denotes the unconditional maximum-likelihood estimate of a parameter, i.e. both the 
parameters of interest and the nuisance parameters are varied to maximise the likelihood function. The 
double circumflex denotes a conditional maximum-likelihood estimate, i.e. an estimate for given fixed 
values of the parameters of interest a. 

Systematic uncertainties and their correlations [107] are modelled by introducing nuisance parameters 
6 described by likelihood functions associated with the estimate of the corresponding effect. System¬ 
atic uncertainties that affect multiple measurements are modelled with common nuisance parameters to 
propagate the effect of these uncertainties coherently to ah measurements. Most experimental systematic 
uncertainties are modelled independently for the yfs = 7 and 8 TeV data samples, reflecting independ¬ 
ent assessments of these uncertainties, but a subset of these uncertainties, e.g. material effects and some 
components of the jet energy scale, are considered common to the two data taking periods and are corres¬ 
pondingly described by a common set of nuisance parameters. 

Components of theoretical uncertainties, scale uncertainties of a given Higgs boson production process 
as well as PDF-induced uncertainties, that affect the inclusive signal rate are described with common 
nuisance parameters in ah channels, whereas components of theoretical uncertainties that affect the ac¬ 
ceptance of individual channels arc modelled with separate nuisance parameters for each decay channel. 
Specifically, since PDF-induced uncertainties and scale uncertainties are described by separate nuisance 
parameters, these uncertainties are effectively treated as uncorrelated. The PDF uncertainties of the in¬ 
clusive rates are treated as correlated for WH , ZH and VBF production, as anti-correlated for gg —» ZH 
and qq —> ZH production and as uncorrelated for ggF and ttH production. A cross check with the full 
correlation matrix as given in Ref. [32] show no differences larger than 1% for the most generic model 
(Section 5.5.3). Similarly, the effects of correlations between Higgs boson branching ratios and partial 
decay widths have been determined to be negligible, and are ignored in the combinations, except for 
the branching ratios to WW* and ZZ* which are treated as fully correlated. When results are provided 
with a breakdown of the systematic uncertainties in experimental and theoretical uncertainties, the the¬ 
oretical uncertainties correspond to the influence of ah nuisance parameters that can affect Higgs boson 
signal distributions, e.g. parton density functions related to Higgs boson production, QCD scale un¬ 
certainties related to Higgs boson production processes and uncertainties on the Higgs boson branching 
ratios. Theoretical uncertainties that exclusively affect background samples are included in the systematic 
uncertainty components. 

The choice of the parameters of interest depends on the test under consideration, with the remaining 
parameters being “profiled", i.e., similarly to nuisance parameters they are set to the values that maximise 
the likelihood function for the given fixed values of the parameters of interest. 

Asymptotically, a test statistic -2 In A(or) of several parameters of interest a is distributed as ax 1 distribu¬ 
tion with n degrees of freedom, where n is the dimensionality of the vector a. In particular, the 100(1 -/?)% 
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confidence level (CL) contours arc defined by -2 In A(ar) < kp, where kp satisfies P(y 2 3 > kp) = />’■ For 
one degree of freedom the 68% and 95% CL intervals arc given by -2 In A(ar) = 1.0 and 4.0, respectively. 
For two degrees of freedom the 68% and 95% CL contours are given by -21nA(ar) = 2.3 and 6.0, re¬ 
spectively. All results presented in the following sections are based on likelihood evaluations and give CL 
intervals under asymptotic approximation. 2 For selected parameters of interest a physical boundary on the 
parameter values is included in the statistical interpretation. For example, branching ratio parameters can 
conceptually not be smaller than zero. The 95% confidence interval quoted for such parameters is then 
based on the profile likelihood ratio restricted to the allowed region of parameter space; the confidence 
interval is defined by the standard x 2 cutoff, which leads to some over-coverage near the boundaries. 

For the measurements in the following sections the compatibility with the Standard Model, psm> is quan¬ 
tified using the p-value J obtained from the profile likelihood ratio Aft* = qtsmX where a is the set of 
parameters of interest and qtsm are their Standard Model values. For a given benchmark coupling model, 
a is the set of Higgs boson coupling scale factors at , and ratios of coupling scale factors Ag probed by that 
model, where the indices i, j refer to the parameters of interest of the model (see Section 5). All other 
parameters are treated as independent nuisance parameters. 


4. Signal-strength measurements 

This section discusses the measurements of the signal-strength parameter g of different production modes 
and decay channels as well as their ratios for a fixed Higgs boson mass hypothesis of mu = 125.36 GeV [23] 
The signal-strength parameter is a measure of potential deviations from the SM prediction under the as¬ 
sumption that the Higgs boson production and decay kinematics do not change appreciably from the SM 
expectations. In particular, the transverse momentum and rapidity distributions of the Higgs boson are as¬ 
sumed to be those predicted for the SM Higgs boson by state-of-the-art event generators and calculations 
of each production process. This assumption is corroborated by studies such as the measurements of dif¬ 
ferential production cross sections [113, 114] and tests of spin and CP properties of the Higgs boson [24, 
115]. 

For the discussion in this section, bbH is assumed to have the same signal strength as ggF, tH the same 
as ttH, and gg —* ZH the same as qq —> ZH, unless noted otherwise. The ggF and bbH processes lead 
to similar event signatures and no attempt is made to separate them in the analyses, thus the assumption 
of equal signal strength implies that the observed ggF signal is interpreted as a mixture of bbH and ggF 
events following their SM ratio of cross sections. The ttH and tH events have similar topologies. The 
gg —» ZH process leads to the same final state as the qq —> ZH process. Whenever WH and ZH are 
combined into VH , their signal strengths arc assumed to be the same. 


4.1. Global signal strength 

In Section 2, the published ATLAS measurements on Higgs boson production and decay modes based on 
individual final states as well as the changes since their publication are summarised. Figure 2 shows the 
updated measurements of the signal-strength parameter g from a simultaneous fit to all decay channels 


2 Whenever probabilities are translated into the number of Gaussian standard deviations the two-sided convention is chosen. 

3 The p-value is defined as the probability to obtain a value of the test statistic that is at least as high as the observed value, 
under the hypothesis that is being tested. 
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analysed, assuming SM values for the cross-section ratios of dilferent Higgs boson production processes 
(or equivalently all /q’s of Eq. (2) are set to be equal). In the fit, the SM predictions of the signal yields are 
scaled by decay-dependent signal-strength parameters, independent of production processes. Compared 
to the separate measurements shown in Fig. 1, small changes are observed, resulting from the assignment 
of the Higgs boson yields in the ttH searches to appropriate decay channels, namely H —> WW*, H —> tt 
and H —> bb.' The central values all increase slightly due to the high observed signal-strength values of 
the ttH searches, but the uncertainties are barely improved because of the limited significance obtained 
for the ttH production process with the current dataset. The most significant change in the signal strength 
is observed for the H —> bb decay. The combination of the VH{—> bb) analysis and the ttH{—> bb) 
search leads to an observed (expected) significance of 1.8 (2.8) standard deviations for the H —> bb decay 
channel. 

Assuming a multiplier common to all decay modes, signal-strength measurements of individual decay 
modes can be combined to give a global and more precise measurement, providing the simplest consist¬ 
ency test with the SM expectation. Combining all measurements using the profile likelihood ratio A (p.) 
results in a global signal-strength value of 

p = 1.18 ^[4 = 1.18 ±0.10 (stat.) ± 0.07 (syst.) (theo.), 

where the labels stat., syst. and theo. refer to statistical, systematic, and signal theoretical uncertainties, 
respectively. The signal theoretical uncertainty includes contributions from uncertainties in SM cross 
sections and branching ratios as well as in the modelling of the production and decays of the Higgs 
boson, as discussed in Section 3. The theoretical uncertainties of background processes are included in 
the uncertainty labelled as systematic uncertainty. 

The uncertainty on the global signal strength has comparable statistical and systematic components and is 
significantly reduced compared to the individual measurements, as illustrated in Fig. 2. Here, the largest 
source of experimental systematic uncertainty is from background estimates in the analyses of individual 
channels. This result is consistent with the SM expectation of p = 1, with a p- value of 18%, All individual 
measurements of the signal-strength parameters are consistent and compatible with the combined value, 
with a p-value of 76%. 

Performing independent combinations of measurements at yfs = 7 and 8 TeV independently lead to 
signal-strength values of 

p(l TeV) = 0.75 + _° 0 H = 0.75 ±°;| (stat.) + _° 0 \l (syst.) + _° 0 ™ (theo.), and 

q(8TeV) = 1.28!°;}^ = 1.28 ±0.11 (stat.)(syst.) (theo.) 

at these two energies. The relative theoretical uncertainty of ~ 7% on the measured p value at -\Js = 8 TeV 
arises predominantly from the uncertainty on the total cross section, but is nevertheless smaller than the 
corresponding uncertainty of ~ 9% on the total SM cross section shown in Table 1, because p is effectively 
a weighted average of the signal-strength measurements in all categories: the contributions from VBF 
and VH production, which have comparatively small theoretical uncertainties, have larger weights in this 
average than in the total cross section. 


4 The measurement of the qqlgg —> ttH signal strength in the multiple-lepton decay mode contributes to all final states with 
leptons in Fig. 2, according to the prediction of MC simulation, i.e. predominantly to the H —> WW* and H —> ZZ* final 
states. 
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Figure 2: The observed signal strengths and uncertainties for different Higgs boson decay channels and their com¬ 
bination for mu = 125.36 GeV. Higgs boson signals corresponding to the same decay channel are combined together 
for all analyses, assuming SM values for the cross-section ratios of different production processes. The best-fit val¬ 
ues are shown by the solid vertical lines. The total ±lcr uncertainties are indicated by green shaded bands, with 
the individual contributions from the statistical uncertainty (top), the total (experimental and theoretical) systematic 
uncertainty (middle), and the signal theoretical uncertainty (bottom) on the signal strength shown as horizontal error 
bars. 


4.2. Individual production processes 

In addition to the signal strengths of different decay channels, the signal strengths of different production 
modes are also determined, exploiting the sensitivity offered by the use of event categories in the analyses 
of all channels. 

The Higgs boson production modes can be probed with four signal-strength parameters: /u gg p, pvm ; , 
and p„//, one for each main production mode, combining Higgs boson signals from different decay 
channels under the assumption of SM values for the ratios of the branching ratios of different Higgs 
boson decays. This assumption is equivalent to set all /u/s in Eq. (2) to be equal. The SM predictions of 
the signal yields are scaled by these four production-dependent parameters. The best-fit values of these 
parameters for the Vs = 8 TeV data separately and in combination with the \[s = 7 TeV data are shown 
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in Table 4. Uncertainty components from statistics, systematics, and signal theory arc also shown. The 
accuracy with which the uncertainties are broken down is limited by the precision of the fit and more 
importantly by the approximations made in individual analyses when neglecting uncertainties which are 
small with respect to, e.g., the statistical uncertainty. The \[s = 7 and 8 TeV combined values with their 
total uncertainties are also illustrated in Fig. 3. The yfs = 7 TeV data are included in the combinations 
only, as they have limited statistical power to distinguish between different production modes. The signal- 
strength measurements are in reasonable agreement with the SM predictions of unity. Although the results 
support the SM prediction of the ttH production (see Section 4.4), this production process remains to be 
firmly established in future LHC runs. Thus, a 95% CL upper limit on its signal strength is also derived. 
Combining the results from various analyses with sensitivity to ttH production, the observed and expected 
limits are [i t m < 3.2 and 1.4, respectively. 


Table 4: Measured signal strengths jj at >hh = 125.36 GeV and their total ±lcr uncertainties for different production 
modes for the ^fs - 8 TeV data and the combination with the \fs = 7 TeV data. The \fs = 7 TeV data do not have 
sufficient statistical power to yield meaningful measurements for individual production modes, but are included in 
the combination. Shown in the square brackets are uncertainty components: statistical (first), systematic (second) 
and signal theoretical (third) uncertainties. These results are derived using SM values for the ratios of branching 
ratios of different Higgs boson decay channels. 


Production 


Signal strength /u at ra# = 125.36 GeV 

process 

V~* = 

8 TeV 

Combined ^fs = 1 and 8 TeV 

ggF 

1 23 +0 - 25 
• -0.21 1 

1 + 0.16 + 0.10 + 0 . 16 ] 
[- 0.16 - 0.08 - 0.11 I 

-t! 

[+ 0.14 + 0.09 + 0 . 16 ] 
[- 0.14 - 0.08 - 0 . 12 ] 


VBF 

1 55 +039 

w - 0.35 | 

[+ 0.32 + 0.17 + 0 . 13 ] 
[- 0.31 - 0.13 - 0.11 | 

1.23 ±0.32 

[+ 0.28 + 0.13 + 0.11 | 
[- 0.27 - 0.12 - 0 . 09 ] 


VH 

0.93 ± 0.39 

[+ 0.37 + 0.20 + 0 . 12 ] 
L — 0.33 - 0.18 - 0 . 06 ] 

0.80 ±0.36 

[+ 0.31 + 0.17 + 0 . 10 ] 
L 0.30 - 0.17 - 0 . 05 ] 


ttH 

1.62 ±0.78 

[+ 0.51 + 0.58 + 0 . 28 ] 
1 - 0.50 - 0.54 - 0.10 1 

1.81 ±0.80 

[+ 0.52 + 0.58 + 0.31 | 
1 - 0.50 - 0.55 - 0.12 1 



Table 5: Measured cross sections of different Higgs boson production processes at = 8 TeV for m H = 
125.36 GeV obtained from the signal-strength values of Table 4. Their SM predictions can be found in Table 1. 
Shown in the square brackets are uncertainty components: statistical (first), systematic (second) and signal theor¬ 
etical (third) uncertainties. The theoretical uncertainties here arise from the modelling of Higgs boson production 
and decays. These results are derived using the SM values of the Higgs boson decay branching ratios. 


Production process Cross section [pb] at yfs = 8 TeV 


ggF 

23.9 ± 3.6 

[+ 3.1 + 1.9 + 1.01 
[- 3.1 - 1.6 - 1 . 0 ] 

VBF 

2.43 ±0.58 

[+ 0.50 + 0.27 + 0 . 19 ] 
[- 0.49 - 0.20 - 0 . 16 ] 

VH 

1.03 ±0.53 

[+ 0.37 + 0.22 + 0 . 13 ] 
L 0.36 - 0.20 - 0 . 06 ] 

ttH 

0.24 ±0.11 

[+ 0.07 + 0.08 + 0.011 
1 - 0.07 - 0.08 - 0.01 1 


The signal-strength measurements shown in Table 4 are extrapolated to total cross-section measurements 
for each production process, as shown in Table 5 for yfs = & TeV, with the further assumption of SM val- 
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Figure 3: The best-fit signal-strength values of different production modes determined from the combined fit to the 
Vs = 7 and 8 TeV data. Higgs boson signals corresponding to the same production process but from different decay 
channels are combined together, assuming SM values for the ratios of the branching ratios of different Higgs boson 
decay channels. The inner and outer error bars correspond to 68% CL and 95% CL intervals. Total uncertainties 
combining statistical, experimental and theoretical systematic uncertainties are shown. 


ues for the Higgs boson decay branching ratios. The theoretical uncertainties on the absolute values of the 
SM Higgs boson production cross sections are thereby removed, but significant theoretical uncertainties 
remain, related to the modelling of the Higgs boson production and of the acceptance of the event selec¬ 
tion. One can sum the different cross sections to obtain an overall extrapolated cross section for Higgs 
boson production. The measurement is performed at sfs -1 TeV as well despite of the limited statistical 
power of the dataset. The resulting total Higgs boson production cross sections at the two energies are 

07/(7 TeV) = 22.1 pb = 22.1 +53 (stat .)+23 (syst.) ^j 4 (theo.) pb, and 

07/(8 TeV) = 27.7 ± 3.7 pb = 27.7 ± 3.0 (stat.) (syst.) (theo.) pb, 

to be compared with the theoretical predictions of 17.4 ± 1.6 pb at yfs = 7 TeV and 22.3 ± 2.0 pb at 
Vx = 8 TeV, as shown in Table 1. 

These cross sections are different from what one would naively expect from the global signal-strength val¬ 
ues discussed in Section 4.1, particularly for \[s = 1 TeV. The differences are largely the result of analysis 
categorisation. Categories often explore production processes or phase-space regions with distinct signal- 
event topologies. The resulting high signal-to-background ratios can significantly improve the precision 
of the signal-strength measurements. However, these categories often account for small fractions of the 
production cross section and thus have limited impact on the total cross-section measurement, which is 
dominated by processes with larger expected cross sections. One good example is the VBF category. It 
contributes significantly to the global signal-strength measurement, but has a relatively minor impact on 
the total cross-section measurement. 
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4.3. Boson and fermion-mediated production processes 

The Higgs boson production processes can be categorised into two groups according to the Higgs boson 
couplings to fermions (ggF and ttH) or vector bosons (VBF and VH). Potential deviations from the SM 
can be tested with two signal-strength parameters, p [ gV+ttH = (// gF = v/,//) and /^bf + vh = (/4bf = 
for each decay channel /, assuming SM values for the ratio of ggF and ttH cross sections and the ratio 
of VBF and VH cross sections. Signal contaminations from one group to another, e.g. ggF events with 
two jets passing the VBF selection, are taken into account in the simultaneous fit. The 68% and 95% CL 
two-dimensional contours of n{ oV+ttH and F vbf+v// °f dlc five main decay channels are shown in Fig. 4. 
The measurements of H —» p/r and H —> Zy decays have relatively poor sensitivities and are therefore 
not included in the figure. The cutoff in the contours of the H-* yy and H —> ZZ* decays is caused by 
the expected sum of signal and background yields in one of the contributing measurements going below 

f 

zero in some regions of the parameter space shown in Fig. 4. The SM expectation of = 1 and 

f 

rvBF+VH = ' ' s within the 68% CL contour of most of these measurements. 




0 




t 

ggF+ttH 


f f 

Figure 4: Likelihood contours in the (p gF+f(W ,Py BF+l/w ) plane for a Higgs boson mass mu = 125.36 GeV measured 
separately for H —> WW*, ZZ*. hb, yy and tt decays. SM values are assumed for the relative contributions between 
ggF and ttH and between VBF and VH production. The straight lower portions of the H —» yy and H —> ZZ* —> 4f 
contours are due to the small numbers of events in these channels and the requirement of a positive probability 
density function. The best-fit values to the data (+) and the 68% (full) and 95% (dashed) CL contours are indicated, 
as well as the SM expectation (★). 


The relative production cross sections of the processes mediated by vector bosons and by fermions can 

f f 

be tested using the ratio r°°F+nH- When measured separately for each decay channel, this ratio 

reduces to the ratio of production cross sections because the Higgs boson decay branching ratios cancel 
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and is equivalent to the ratio of /q defined in Section 4.1, i.e., 


^VBF +VH _ °~VBF +VHICgg¥+ttH _ RWBF+VH _ ^ ^ 

V f ggF+t,H [ <T VBF +VH/CT gg F +tt H \ SM VggF+ttH 

The observed ratios are shown in Table 6 and illustrated in Fig. 5 for the five main decay channels. The 

f 

signal-strength parameter n gg¥+ttH of each decay channel is profiled in the fit. The combination of these 
measurements yields an overall value of the ratio of cross sections for the vector-boson- and fermion- 
mediated processes (relative to its SM prediction): 

^Combined = 0.96 ^ = 0.96 (stat.) ± 0 J 3 (syst.) (theo.). 


Table 6: The best-fit values and their uncertainties for the ratio Rff of cross sections for the vector-boson- and 
fermion-mediated production processes relative to their SM values at mu = 125.36 GeV for the individual decay 
channels and their combination. Shown in the square brackets are uncertainty components: statistical (first), sys¬ 
tematic (second) and signal theoretical (third) uncertainties. These results are independent of the Higgs boson decay 
branching ratios. 


Decay channel 
H —* yy 
H —> ZZ* 

H -> WW* 

H —» TT 
H —> bb 
Combined 


Cross-section ratio Rff 


0.56 

0.18 

1.47 

0.81 

0.33 

0.96 


+ 0.66 1 

[+ 0.62 

+ 0.15 

+ 0.18 

- 0.45 | 

[- 0.42 

- 0.09 

- 0.15 

+ 1.20 1 

[+ 1.16 

+ 0.23 

+ 0.23 

- 0.52 | 

[- 0.50 

- 0.05 

- 0.16 

+ 0.80 

[+ 0.63 

+ 0.37 

+ 0.31 

- 0.54 | 

[- 0.47 

- 0.19 

- 0.18 

+ 2.19 1 

[+ 1.36 

+ 1.68 

+ 0.39 

- 0.49 | 

[- 0.41 

- 0.15 

- 0.23 

+ 1.03 1 

[+ 0.39 

+ 0.94 

+ 0.18 

- 0.25 1 

1 - 0.20 

- 0.14 

- 0.06 


+ 0.43 + 0.33 + 0.20 + 0.18 
- 0.31 - 0.26 - 0.13 - 0.10 
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Figure 5: The ratios of cross sections for the vector-boson- and fermion-mediated processes relative to their SM 
values at % = 125.36 GeV, measured in the individual Higgs boson decay final states and their combination, 
/'’combined (see text). The inner and outer error bars represent 68% CL and 95% CL intervals, combining statistical 
and systematic uncertainties. These measurements are independent of Higgs boson decay branching ratios. 
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4.4. Ratios of production cross sections and partial decay widths 


At the LHC, the Higgs boson production cross sections and decay branching ratios cannot be separately 
determined in a model-independent way as only their products are measured. However, the ratios of 
cross sections and ratios of branching ratios can be disentangled without any assumptions, within the 
validity of the narrow width approximation of the Higgs boson. By normalising to the cross section of the 
gg —> H —» WW* production process, cr(gg —>//—> WW*), the yields of other Higgs boson production 
modes and decay channels can be parameterised using the ratios of cross sections and ratios of branching 
ratios. For the production and decay i —* H —* f, the yield is then 


CTj ■ 



= <r(gg -> H 


WW*)x 



x 


_£/_ 

r ww 


( 6 ) 


The ratio of branching ratios in the above equation is substituted by the equivalent ratio of partial decay 
widths. The ratios extracted from the measured yields are independent of theoretical predictions on the 
inclusive cross sections and partial decay widths (and thus branching ratios). Furthermore, many exper¬ 
imental systematic uncertainties cancel in the ratios. The residual theoretical uncertainties are related to 
the modelling of the Higgs boson production and decay, which impacts the signal acceptance calculations. 
The gg —» H —> WW* process is chosen as the reference because it has both the smallest statistical and 
overall uncertainties, as shown in Fig. 2. 

The yfs = 1 and 8 TeV data are fitted with <r(gg —> H —> WW*), cr,-/cr gg p and Tf/F^w* as parameters 
of interest and the results are listed in Table 7, together with the SM predictions [32]. The results after 
normalising to their SM values are illustrated in Fig. 6. The results of <r(gg —> H —> WW*) and o- ; /cr gg F 
from the combined analysis of the \[s ~ 1 and 8 TeV data are shown for yfs = 8 TeV, assuming the SM 
values for cr,(7 TeV)/cr,(8 TeV). The WH and ZH production processes are treated independently in the 
fit to allow for direct comparisons with theoretical predictions. The searches for H —» pp and H —» Zy 
decays are included in the fit, but the current datasets do not result in sensitive measurements for these 
two decays. Therefore only 95% CL upper limits are derived, namely 0.006 for F^/Tww and 0.078 
for Yzy/Tww- The p-value of the compatibility between the data and the SM predictions is found to be 
80%. 

The results exhibit a few interesting features that are worth mentioning. As a multiplicative factor com¬ 
mon to all rates in this parameterisation, cr(gg —>//—> WW*) is pulled up in the fit to accommodate the 
observed large global signal-strength value (Section 4.1). The best-fit value of cr(gg —>//—> WW*) is 
approximately 15% above the SM prediction, to be compared to the significantly lower value of 0.98 
found from the stand-alone measurement from the H —» WW* decay (see Fig. 1). Moreover, there are by 
construction large anti-correlations between cr(gg —>//—> WW*), <x,V<x gg F and Vj/V ww- 

Table 8 shows the observed and expected significances in units of standard deviations of the VBF, WH, 
ZH and ttH production processes. Listed under VH are the combined significances of WH and ZH 
production, assuming the SM value for their relative cross sections. The significance is calculated from 
a likelihood scan, where the contributions from other processes are fixed at their best-fit values. As the 
gg —» H —> WW* process is chosen as the reference, the significances are calculated using the observable 
cr(gg —> H —> WW" ) for the ggF process and the cross-section ratios c,7cr ggF for all other processes. The 
cross-section ratios are independent of the Higgs boson decay branching ratios and have the advantage of 
the cancellation of many experimental uncertainties. The result provides an unequivocal confirmation of 
the gluon fusion production of the Higgs boson with its significance exceeding well above five standard 
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Table 7: Best-fit values of cr(gg —>//—> WW*), cr,/<x ggF and T//T ww . for a Higgs boson with m H — 125.36 GeV 
from the combined analysis of the ^ = 1 and 8 TeV data. The cross-section ratios are given for yfs = 8 TeV as¬ 
suming the SM values for <x,(7 TeV)/cr,-(8 TeV). Shown in square brackets are uncertainty components: statistical 
(first), systematic (second) and signal theoretical (third) uncertainties. The SM predictions [32] are shown in the 
last column. 


Parameter Best-fit value SM prediction 


-> H -> WW*) (pb) 

4.86 | 

[+0.76 +0.52 +0.26] 
[-0.74 -0.48 -0.18] 

4.22 ± 0.47 

O' VBF/cr ggF 

0-081 tS | 

[+0.031 +0.016 +0.008] 
[-0.024 -0.010 -0.005] 

0.082 ± 0.009 

O' WH/cr gg F 

0053 to* 

[+0.032 +0.018 +0.008] 
L-0.023 -0.012 -0.004] 

0.036 ± 0.004 

CZff/cTggF 

Ooi3+ 0 - 030 
u.uij _ () 014 I 

[+0.021 +0.020 +0.005] 
[-0.013 -0.005 -0.002] 

0.021 ± 0.002 

ctm/c ggF 

0 012+ 0 - 007 
-0.005 

[+0.005 +0.004 +0.0014] q AQ-y + Q Q01 

[-0.004-0.003 -0.0005] u.uu/ + u.uui 

r rr /r ww . 

ooio +0 - 003 

u.uiu _ () 0 03 

[+0.003 +0.002 +0.0006 ] 0 01036 + 0 00011 

[-0.002 -0.001 -0.0004] U ' U “ 1 

r zz ./r ww 

0 15 +0 - 05 

U - 1J -0.04 

[+0.046 +0.022 +0.008] 
[-0.036 -0.013 -0.005] 

0.124 ±<0.001 


p 

1 + 
p p 

[+0.112 +0.084 +0.032] 
[-0.090 -0.053 -0.017] 

0.285 ± 0.006 

T'bbl^'wW 

1 53 +L64 

1 -0.94 [ 

[+1.17 +1.11 +0.30] 
1-0.69 -0.63 -0.12 | 

2.60 ±0.12 


Table 8: The observed and expected significances in units of standard deviations for different Higgs boson produc¬ 
tion processes except ggF production which is well established (see text). The significances of VH production are 
obtained by combining the WH and ZH processes, assuming the SM value for their relative cross sections. All 
significances are calculated under the asymptotic approximation [112], 


Process 

VBF 

ttH 

WH 

ZH 

VH 

Observed 

4.3 

2.5 

2.1 

0.9 

2.6 

Expected 

3.8 

1.5 

2.0 

2.1 

3.1 
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Figure 6: The gg —> H —> WW* cross section, ratios of cross sections and of partial decay widths relative to their 
SM values at mg = 125.36 GeV from the combined analyses of the \fs = 7 and 8 TeV data. The inner and outer 
error bars on the measurements are 68% CL and 95% CL intervals. The SM predictions are shown as the vertical 
line at unity with grey bands representing theoretical uncertainties on the ratios of inclusive cross sections and of 
partial decay widths. 


deviations. Furthermore, the result also offers strong evidence, at 4.3 standard deviations, of vector-boson 
fusion production and supports the SM assumptions of production in association with vector bosons or a 
pair of top quarks. 

An alternative parameterisation normalising the ratios of cross sections and of branching ratios to their 
SM values is presented in Appendix A. 
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5. Coupling-strength fits 


f 

In the previous section signal-strength parameter p for a given Higgs boson production or decay mode 
is discussed. For a measurement of Higgs boson coupling strengths, production and decay modes cannot 
be treated independently, as each observed process involves at least two Higgs boson coupling strengths. 
Scenarios with a consistent treatment of coupling strengths in production and decay modes are studied in 
this section. All uncertainties on the best-fit values shown take into account both the experimental and the¬ 
oretical systematic uncertainties. For selected benchmark models a breakdown of parameter uncertainties 
in statistical uncertainties and in experimental and theoretical systematic uncertainties is presented. 


5.1. Framework for coupling-strength measurements 

Following the leading-order (LO) tree-level-motivated framework and benchmark models recommended 
in Ref. [32], measurements of Higgs boson coupling-strength scale factors kj arc implemented for the 
combination of all analyses and channels summarised in Table 3. 


5.1.1. Structure and assumptions of the framework for benchmark models 


The framework is based on the assumption that the signals observed in the different channels originate 
from a single narrow resonance with a mass near 125.36 GeV. The case of several, possibly overlapping, 
resonances in this mass region is not considered. Unless otherwise noted, the Higgs boson production 
and decay kinematics are assumed to be compatible with those expected for a SM Higgs boson, similar 
to what was assumed for the signal-strength measurements of Section 4. 

The width of the assumed Higgs boson near 125.36 GeV is neglected in the Higgs boson propagator, i.e. 
the zero-width approximation is used. In this approximation, the cross section —> H —» /) for on-shell 

measurements can always be decomposed as follows: 


cr(i 


H 


/) = 


Ci(Kj) ■ r f(Kj) 

r H {Kj) 


(V) 


where o"j is the Higgs boson production cross section through the initial state i, Vj its the partial decay 
width into the final state / and T# the total width of the Higgs boson. The index j runs over all Higgs 
boson couplings. The components of a, Tf, and Vn of Eq. (7) are expressed in scale factors kj of the 
Higgs boson coupling strengths to other particles j that are motivated by the leading-order processes that 
contribute to production or decay, and are detailed in Section 5.1.2. All scale factors are defined such 
that a value of Kj = 1 corresponds to the best available SM prediction, including higher-order QCD and 
EW corrections. This higher-order accuracy is generally lost for kj + 1, nevertheless higher-order QCD 
corrections approximately factorise with respect to coupling rescaling and are accounted for wherever 
possible. 

Modifications of the coupling scale factors change the Higgs boson width T h{Kj) by a factor k 2 h {kj) with 
respect to the SM Higgs boson T™, 


r H (Kj) = k 2 h (kj ) • r™, 
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where k 2 h {kj) is the sum of the scale factors k 2 weighted by the corresponding SM branching ratios. The 
total width of the Higgs boson increases beyond modifications of kj if invisible or undetected Higgs boson 
decays 5 occur that are not present in the SM. Including a Higgs boson branching fraction BRi u to such 
invisible or undetected decays, the full expression for the assumed Higgs boson width becomes 


!*(*,•, BRl, u .) 


K 2 H (Kj) 

(l-BRi., u .) 


r SM 

1 H ' 


( 8 ) 


As BRi u scales all observed cross-sections of on-shell Higgs boson production cr(i —» H —> /), some 
assumption about invisible decays must be made to be able to interpret these measurements in terms of 
absolute coupling-strength scale factors Kj. The signal-strength measurements of off-shell Higgs boson 
production [21], on the other hand, is assumed to only depend on the coupling-strength scale factors and 
not on the total width [103, 104], i.e. 


cr off (i -» ~ /cjoff ' ( 9 ) 

where the additional assumption of non-running coupling-strength scale factors, /r / 0 ir = K J(m allows r// to 
be constrained using using Eq. (8), from a simultaneous measurement of on-shell and off-shell measure¬ 
ments. While this assumption of non-running coupling-strength scale factors cannot hold universally for 
ggF and VBF production without violating unitarity, it is assumed to hold in the region of phase space of 
the off-shell H* —» WW and H* —> ZZ measurements described in Section 2.9 which is relatively close to 
the on-shell regime [116]. Alternatively, ratios of coupling-strength scale factors can be measured without 
assumptions on the Higgs boson total width, as the identical contributions of T# to each coupling strength 
cancel in any ratio of these. 

Finally, only modifications of coupling strengths, i.e. of absolute values of coupling strengths, are taken 
into account, while the tensor structure of the couplings is assumed to be the same as in the SM. This 
means in particular that the observed state is assumed to be a CP-even scalar as in the SM. This assumption 
was tested by both the ATFAS [24] and CMS [115] Collaborations. 


5.1.2. Characterisation of the input measurements in terms of coupling strengths 

The combined input channels described in Table 3 probe eight different production processes: cr(ggF), 
cr(VBF), cr(WH), cr(qq —> ZH ), <r(gg —> ZH), cr(bbH), cr{ttH ), and cr(tH) whose SM cross sections are 
listed in Table l. 6 Table 9 summarises the Higgs boson coupling-strength characteristics of all production 
processes and lists the rate scaling behaviour in terms of Higgs boson coupling-strength scale factors. 

The ggF production process (Fig. 7a) involves a loop process at lowest order, with contributions from 
t- and 6-quark loops and a small interference between them. The VBF production (Fig. 7b) process 
probes a combination of Kw and Ky coupling-strength scale factors, with a negligible amount («: 0.1%) 
of interference between these tree-level contributions. 

The qq —> WH and qq —> ZH processes (Fig. 8a) each probe a single coupling strength, with scale 
factors kw and kz, respectively. The gluon-initiated associated production of a Higgs boson with a Z 
boson, cr(gg —> ZH), is characterised by gluon-fusion-style production involving t, 6-quark loops where 

5 Invisible final states can be directly searched for through the £“ lss signature [28], An example of an undetected mode would 
be a decay mode to multiple light jets, which presently cannot be distinguished from multijet backgrounds. 

6 The ZH production cross section quoted in Table 1 comprises both the qq —> ZH and gg —> ZH processes. 
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Figure 7: Feynman diagrams of Higgs boson production via (a) the ggF and (b) VBF production processes. 



Figure 8: Feynman diagrams of Higgs boson production via (a) the qq —* VH and (b,c) gg —> ZH production 
processes. 


the Z boson is always radiated from the fermion loop and the Fliggs boson is either radiated directly 
from the fermion loop (Fig. 8b), or is radiated from the outgoing Z boson (Fig. 8c). The cross section of 
gg — * ZH production is sensitive to the relative sign between K t and kz due to interference between these 
contributions. This separate treatment of gg —> ZH production is not present in the framework described 
in Ref. [32], 

The ttH production process (Fig. 9a) directly probes the Higgs boson coupling strength to top quarks, 
parameterised in the framework with the scale factor K t . Tree-level tH production, comprising the pro¬ 
cesses qg —* tHbq' (Fig. 9b, 9c) and gb —> WtH (Fig. 9d, 9e), is included as background to events in all 
reconstructed ttH categories, and has for SM Higgs boson coupling strengths a large destructive interfer¬ 
ence [69] between contributions where the Higgs boson is radiated from the W boson and from the top 
quark. The SM cross section for tH production is consequently small, about 14% of the ttH cross section. 
However, for negative K t the interference becomes constructive and, following Table 9, the cross section 
increases by a factor of 6 (13) for | K t \ - \ Kw I = 1 for the gb —> WtH {qg —> tHbq') process, making the 
tH process sensitive to the relative sign of the W and top-quark coupling strength, despite its small SM 
cross section. The modelling of tH production is not present in the framework described in Ref. [32], 

The bbH (Fig. 9a) production process directly probes the Higgs boson coupling strength to 6-quarks, with 
scale factor ki,. Simulation studies using bbH samples produced in the four-flavour scheme [82, 96] have 
shown that the ggF samples are a good approximation for bbH production for the most important analysis 
categories, therefore bbH production is always modelled using simulated ggF events (see Section 2.10). 

The combined input channels probe seven Higgs boson decay modes. Five of these decay modes, H —> 
WW*, H —> ZZ*, H —> bb, H —» rr, and H —* g/i each probe a single coupling-strength scale factor to 
either a gauge boson (Fig. 10a) or to a fermion (Fig. 10b). The remaining two decay modes, H —> yy and 
H —> Zy are characterised by the interference between W boson or top-quark loop diagrams (Fig. 11). 
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Figure 10: Feynman diagrams of Higgs boson decays (a) to W and Z bosons and (b) to fermions. 


These modes probe the W and t coupling strengths as well as their relative sign through interference 
effects. 



(a) (b) (c) 

Figure 11: Feynman diagrams of Higgs boson decays to a pair of photons, or to a photon and a Z boson. 


For completeness it should be noted also that the ggF, tH and gg —> ZH cross sections expressed in Fliggs 
boson coupling strengths depend on the kinematic selection criteria used. The b-t interference expression 
quoted in Table 9 for ggF is valid for the inclusive cross section, but in events with additional jets the 
top-quark loop dominates, and the observed interference is somewhat smaller. For gg —> ZH production 
the effect of phase-space dependence was estimated for H —> bb decays with a variant of the coupling 
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model that introduces separate coupling-dependent cross-section expressions for each of the Z boson 
p-\- bins of the H —> bb analysis. The effect on coupling strength measurements of approximating the 
gg —> ZH production cross section with an inclusive expression instead of using the set of px-dependent 
expressions was determined to be negligible at the current experimental precision, with the largest effect 
being a ~ 0.1 cr reduction of the expected sensitivity in the determination of the relative sign of the W/Z 
couplings. Neither this phase-space dependence, nor that of ggF are considered in this paper. For the 
tH process on the other hand, which features a comparatively large W-t interference term, the effect of 
phase-space dependence is taken into account, even though Table 9 only lists the inclusive expression. 


5.1.3. Effective coupling-strength scale factors 

In some of the fits, effective scale factors K g , K y and Ky y are introduced to describe the processes gg —* H, 
H —> yy and H —> Zy, which are loop-induced in the SM, as shown in Figures 7a and 11, respectively. 
In other fits they are treated as a function of the more fundamental coupling-strength scale factors K t , Kb, 
kw, and similarly for all other particles that contribute to these SM loop processes. In these cases, the 
loop contributions arc expressed in terms of the fundamental coupling strengths, including all interference 
effects, as listed for the SM in Table 9. The loop process gg —> ZH is never treated as an effective scale 
factor, as unlike in the other loop processes, a ggHZ contact interaction from new physics would likely 
show a kinematic structure very different from the SM gg —» ZH process [58] assumed in the current 
study and is expected to be suppressed. What then remains of BSM effects on the gg —> ZH process are 
modifications of the Fliggs boson couplings to the top quark (Fig. 8b) and the Z boson (Fig. 8c), which 
are taken into account within the limitation of the framework by the coupling-strength scale factors k, and 
«z- 

5.1.4. Strategies for measurements of absolute coupling strengths 

As all observed Fliggs boson cross sections in the LO framework are inversely proportional to the Fliggs 
boson width (Eq. (7)), which is not experimentally constrained to a meaningful precision at the LHC, 
only ratios of coupling strengths can be measured at the LHC without assumptions about the Higgs boson 
width. To make measurements of absolute coupling strengths, an assumption about the Higgs boson 
width must be introduced. 

The simplest assumption is that there are no invisible or undetected Higgs boson decays, i.e. BRj u = 0 is 
assumed in Eq. (8). An alternative, less strong assumption, is that Kw < 1 and Ky < I [32]. This assumption 
is theoretically motivated by the premise that the Higgs boson should solve the unitarity problem in vector 
boson scattering and also holds in a wide class of BSM models. In particular, it is valid in any model 
with an arbitrary number of Higgs doublets, with and without additional Higgs singlets. The assumption 
is also justified in certain classes of composite Higgs boson models. A second alternative is to assume 
that the coupling strengths in off-shell Higgs boson production are identical to those for on-shell Higgs 
boson production. Under the assumption that the off-shell signal strength and coupling-strength scale 
factors are independent of the energy scale of Higgs boson production, the total Higgs boson decay width 
can be determined from the ratio of off-shell to on-shell signal strengths [21]. The constraint BR L U > 0, 
motivated by the basic assumption that the total width of the Higgs boson must be greater or equal to 
the sum of the measured partial widths, always introduces a lower bound on the Higgs boson width. The 
difference in effect of these assumptions is therefore mostly in the resulting upper limit on the Higgs boson 
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Table 9: Overview of Higgs boson production cross sections cr,, the Higgs boson partial decay widths T/ and 
the Higgs boson total width T//. For each production or decay mode the scaling of the corresponding rate in 
terms of Higgs boson coupling-strength scale factors is given. For processes where multiple amplitudes contribute, 
the rate may depend on multiple Higgs boson coupling-strength scale factors, and interference terms may give 
rise to scalar product terms KiK j that allow the relative sign of the coupling-strength scale factors k ,■ and kj to be 
determined. Expressions originate from Ref. [32], except for cr{gg —* ZH ) (from Ref. [58]) and cr(gb — > Will) and 
cr{qb —> tllq' ) (calculated using Ref. [82]). The expressions are given for a/v = 8 TeV and in H = 125.36 GeV and 
are similar for y[s = 1 TeV. Interference contributions with negligible magnitudes have been omitted in this table. 


Production 

Loops 

Interference 

Expression in fundamental coupling-strength scale factors 

otggF) 

/ 

b-t 

K 2 ~ 

K g 

1.06 • Kf + 0.01 • Kfr - 0.07 • K t Kb 

cr(VBF) 

- 

- 

~ 

0.74 • K 2 W + 0.26 • 4 

cr(WH) 

- 

- 

~ 

K 2 

K w 

cr(qq —> ZH) 

- 

- 

~ 

K 2 

K z 

cr(gg -> ZH) 

/ 

Z-t 

K 2 

ggZH 

2.27 -4 +0.37 -k 2 - 1.64 ■ k z k, 

cr(bbH) 

- 

- 

4 

cr{ttH) 

- 

- 

~ 

4 

cr(gb —> WtH) 

- 

W-t 

~ 

1.84 ■ k 2 + 1.57 ■ 4^ — 2.41 • k,kw 

cr(qb —> tHq') 

- 

W-t 

~ 

3.4 • 4 + 3.56 • 4/ 5.96 • k,kw 

Partial decay width 

r bb 

- 

- 

~ 

4 

r ww 

- 

- 

~ 

% 

Fzz 

- 

- 

~ 

4 

Trr 

r 

1 m 

Tyy 

/ 

W-t 

K y ~ 

4 

K 2 

K n 

1.59 • 4 , + 0.07 • 4 - 0.66 • kwKi 

T Z y 

/ 

W-t 

K Zy ~ 

1.12 ■ 4, + 0.00035 • 4 - 0.12 • KwKt 

Total decay width 

r H 

/ 

W-t 

b-t 

K 2 ~ 

k h 

0.57 • 4 + 0.22 • 4 , + 0.09 • k 2 + 

0.06 • 4 + 0.03 • 4 + 0-03 • 4 + 

0.0023 • 4 + 0.0016 • 4 y + 0.00022 • k 2 


width. The assumptions made for the various measurements are summarised in Table 10 and discussed in 
the next sections together with the results. 

5.2. Fermion versus vector (gauge) coupling strengths 

Benchmark coupling models in this section allow for different Higgs boson coupling strengths to fermions 
and bosons, reflecting the different structure of the interactions of the SM Higgs sector with gauge bosons 
and fermions. It is always assumed that only SM particles contribute to the gg —> H, H—> yy, H —> Zy 
and gg —> ZH vertex loops, and modifications of the coupling-strength scale factors for fermions and 
vector bosons are propagated through the loop calculations. Models with and without assumptions about 
the total width are presented. 
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5.2.1. Assuming only SM contributions to the total width 


In the first benchmark model no undetected or invisible Higgs boson decays are assumed to exist, i.e. 
BRi u = 0. The universal coupling-strength scale factors Kp for all fermions and Ky for all vector bosons 
are defined in this model as: 


Ky = Kw = Kz 

Kp = Kt - Kb - Kt - K g - K m . 

As only SM particles are assumed to contribute to the gg —> H loop in this benchmark model, the gluon 
fusion process depends directly on the fermion scale factor k 2 f . Only the relative sign between Kp and Ky 
is physical and hence in the following only Ky > 0 is considered, without loss of generality. Sensitivity to 
this relative sign is gained from the negative interference between the loop contributions of the W boson 
and the /-quark in H—> yy and H —» Zy decays and in gg —> ZH production, as well as from the tH 
processes (see the corresponding expressions in Table 9). 

Figure 12 shows the results of the fits for this benchmark model. Figure 12a illustrates how the decays 
H—> yy, H —> ZZ*, H —> WW*, H —» rr and H —* bh contribute to the combined measurement. The 
slight asymmetry in Kp for H —> WW* and H —* bh decays is introduced by the small contributions of the 
tH and gg —> ZH production processes that contribute to these decay modes, and which are sensitive to 
the sign of Kp due to interference effects. The strong constraint on Kp from H —» WW* decays is related 
to the 3.2cr observation of the VBF production process in this channel[ll]. Outside the range shown in 
Fig. 12a there are two additional minima for H—> yy. The long tails in the H —> bh contour towards high 
values of \Ky\ are the result of an asymptotically disappearing sensitivity of the observed signal strength in 
the bb final states to Ky at large values of Ky. The combined measurement without overlays is also shown 
in Figure 12b. 

Figures 12a and 12b only show the SM-like minimum with a positive relative sign, as the local minimum 
with negative relative sign is disfavoured at the 4.0cr level, which can been seen in the wider scan of 
Kp, where Ky is profiled, shown in Fig. 12c. The likelihood as a function of Ky, profiling Kp, is given in 
Fig. 12d. Around Ky = 0.8 the sign of the chosen profiled solution for Kp changes, causing a kink in the 
likelihood. The profile likelihood curves restricting Kp to either positive or negative values are also shown 
in Fig. 12d as thin curves, and illustrate that this sign change in the unrestricted profile likelihood is the 
origin of the kink. 

Both Kp and Ky are measured to be compatible with their SM expectation and the two-dimensional com¬ 
patibility of the SM hypothesis with the best-fit point is 41%. The best-fit values and uncertainties are: 

Ky = 1.09 ± 0.07 [!ao5(stat-)!S:o3(syst-)!ao3(theo.)] 

k f = 1.11 ±0.16 [!ali(stat.) ^(syst-) !ao5(theo.)]. 


5.2.2. Allowing for invisible or undetected Higgs boson decays in the total width 

The second benchmark model of this section allows for the presence of invisible or undetected Higgs 
boson decays by introducing BR L lL as a free parameter in the expression of Eq. (8) for the Higgs boson 
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(b) (c) (d) 

Figure 12: Results of fits for the two-parameter benchmark model defined in Section 5.2.1 that probes different 
coupling-strength scale factors for fermions and vector bosons, assuming only SM contributions to the total width: 
(a) results of the two-dimensional fit to Kf and Ky, including 68% and 95% CL contours; overlaying the 68% 
CL contours derived from the individual channels and their combination; (b) the same measurement, without the 
overlays of the individual channels; (c) the profile likelihood ratio as a function of the coupling-strength scale 
factors Kp (Ky is profiled) and fd) as a function of k v ( kf is profiled). The dashed curves in (c) and fd) show the 
SM expectations. In (d) the sign of the chosen profiled solution for k f changes at k v ~ 0.8 , causing a kink in the 
likelihood. The profile likelihood curves restricting k f to be either positive or negative are also shown to illustrate 
that this sign change in the unrestricted profile likelihood is the origin of the kink. The red (green) horizontal 
line indicates the value of the profile likelihood ratio corresponding to a 68% (95%) confidence interval for the 
parameter of interest, assuming the asymptotic ;p 2 distribution for the test statistic. 
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total width. The free parameters of this model thus are Kp, Ky and BR I U . Loop processes are still assumed 
to have only SM content. 

With the introduction of BRj u as a free parameter, the assumed Higgs boson width has no intrinsic upper 
bound and an additional constraint must be imposed on the model that infers an upper bound on T#. Both 
choices of constraints on the total width discussed in Section 5.1 are studied: Ky < 1 and v on = K 0 fi. 


(95% CL) K V > 0.93 

K V = 1.13lo|? 
Ky = 1.09 ± 0.07 

k f = 1.05± 0.16 

K . _i -i 7+0.25 

k F~ 1 • 1 '- 0.16 

k f = 1.1 1 ±0.16 

(95% CL) BR IU <0.13 

(95% CL) BRi'U. <0.52 


-sm - "I -07_ 0 21 

H I ~H Qo+i -35 

= I ■O0_o.3i 

■* qq + 0.30 
jm = I •^-'J-0.26 


ATLAS 


y/s = 7TeV, 4.5 - 4.7fb -1 
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68% CL: 
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K V < 1 K on = K off BR i u = 0 
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m H = 125.36 GeV 
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Parameter value 


Figure 13: Results of fits for benchmark models that probe for potential extra contributions to the total width, but 
do not allow contributions from non-SM particles in the H— » yy, gg —> H and H —> Zy loops, with free gauge and 
fermion coupling-strength scale factors Ky, Kp- The estimated values of each parameter under the constraint k v < 1, 
K 0 n = v () ir or BR I U = 0 are shown with markers in the shape of box, circle, or diamond, respectively. The inner and 
outer bars correspond to 68% CL and 95% CL intervals. The confidence intervals of BRj u and, in the benchmark 
model with the constraint Ky < 1, also Ky, are estimated with respect to their physical bounds, as described in the 
text. The numerical values of the fit under the constraint Ky < 1 are shown on the left. Values for the two alternative 
constraints are also shown (in a reduced font size due to space constraints). 

Figure 13 shows the results of fits for this benchmark scenario. For comparison the results of the bench¬ 
mark model of Section 5.2.1 are included, corresponding to the condition BR L LI = 0. The coupling- 
strength scale factors Kp and Ky are measured to be compatible with the SM values and a limit is set on 
the fraction of Higgs boson decays to invisible or undetected final states. The three-dimensional compat¬ 
ibility of the SM hypothesis with the best-fit point is 99% (29%), when applying the Ky < 1 (off-shell) 
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constraint, respectively. When imposing the physical constraint BRj u > 0, the 95% CL upper limit is 
BRj u < 0.13 (0.52), when applying the constraint Ky < 1 (/r on = A: 0 (f). The corresponding expected limit 
on BR, u , under the hypothesis of the SM, is 0.24 (0.71). 

Also shown in Fig. 13 is the uncertainty on the total width that the model variants allow, expressed as 
the ratio T#/r^ M . These estimates for the width are obtained from alternative parameterisations of these 
benchmark models, where the coupling-strength scale factor Kp is replaced by the expression that results 
from solving Eq. (8) for Kp, introducing T///T^ M as a parameter of the model. Figure 13 shows that the 
upper bound on the Fliggs boson width from the assumption /r () n = K (m is substantially weaker than the 
bound from the assumption Ky < 1. These choices of constraints on the Higgs boson width complement 
each other in terms of explored parameter space: the present limit of p 0 s < 5.1 [21] in the combined 
off-shell measurement in the H — > WW* and H —> ZZ* channels effectively constrains Ky to be greater 
than one in the combined tit when exploiting the assumption /c on = K a ff. 

The parameterisation of the off-shell signal strength q () ir in terms of couplings implicitly requires that 
q 0 ff > 0 (see Ref. [21] for details). This boundary condition causes the distribution of the test statistic 
to deviate from its asymptotic form for low values of cr 0 ff, with deviations in p- values of up to 10% 
for cr 0 ff ~ 2.5, which corresponds to the value of cr 0 ff at the upper boundary of the 68% asymptotic 
confidence interval of T#/T™. The upper bound of the 68% CL interval for the scenario /c 0 ff = /c on 
shown in Fig. 13 should therefore be considered to be only approximate. Since the lower bound on 
r H /T s ™ is always dominated by the constraint BR 1U > 0, it is not affected by this deviation from the 
asymptotic behaviour. 


5.2.3. No assumption about the total width 

In the last benchmark model of this section no assumption about the total width is made. In this model 
only ratios of coupling-strength scale factors are measured, choosing as free parameters 

Apy = Kp/Ky 
Kyy - Ky ■ Ky/Kn, 

where Apy is the ratio of the fermion and vector boson coupling-strength scale factors, Kyy is an overall 
scale that includes the total width and applies to all rates, and kh is defined in Table 9. 

Figure 14 shows the results of this fit. Both ratio parameters are found to be consistent with the SM 
expectation and the two-dimensional compatibility of the SM hypothesis with the best-fit point is 41%. 
The best-fit values and uncertainties, when profiling the other parameter, are: 

Ap V = 1.02^;}^[^;}}(stat.)!^(syst.)^(theo.)] 

Kyy = 1.07!g;}^[!g;J{(stat.)!g;^(syst.)!g;^(lheo.)]. 

Similar to the model described in Section 5.2.1, Fig. 14a shows the determination of the sign of Apy 
disfavouring Apy = -1 at approximately 4.0cr, while Fig. 14b shows the two-dimensional likelihood 
contour. The estimates of the two parameters are anticorrelated because only their product appears in the 
model. 
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Figure 14: Results of fits for the two-parameter benchmark model defined in Section 5.2.3 that probes different 
coupling-strength scale factors for fermions and vector bosons without assumptions about the total width: (a) pro¬ 
file likelihood ratio as a function of the coupling-strength scale factor ratio Apy ( Kyy is profiled). The dashed curve 
shows the SM expectation. The red (green) horizontal line indicates the value of the profile likelihood ratio corres¬ 
ponding to a 68% (95%) confidence interval for the parameter of interest, assuming the asymptotic x 2 distribution 
for the test statistic, (b) Results of the two-dimensional fit to Kyy and Apy, including 68% and 95% CL contours. 


5.3. Probing relations within the fermion coupling sector 

The previous sections assumed universal coupling-strength scale factors for all fermions, while many ex¬ 
tensions of the SM predict deviations from universality within the fermion sector [32]. In this section, 
benchmark models are explored that probe the relations between the up- and down-type fermions and 
between the lepton and quark sectors, using the information in the currently accessible channels, in par¬ 
ticular in H —* bb, H —> rr and H —> p// decays and ttH production. The models considered assume that 
only SM particles contribute to the gg —> H, H-* yy, H —> Zy and gg —» ZH vertex loops, and modific¬ 
ations of the coupling-strength scale factors are propagated through the loop calculations. As only ratios 
of coupling-strength scale factors are explored, no assumptions on the total width are made. 


5.3.1. Probing the up- and down-type fermion symmetry 

Many extensions of the SM contain different coupling strengths of the Higgs boson to up-type and down- 
type fermions. This is for instance the case for certain Two-Higgs-Doublet Models (2HDM) [117-119]. 
In this benchmark model the ratio Aj u of down- and up-type fermions coupling-strength scale factors is 
probed, while vector boson coupling-strength scale factors are assumed to be unified and equal to Ky. The 
indices u, d stand for all up- and down-type fermions, respectively. The free parameters are: 

drfu = KrflKu 

Avu = k v /k u 
K m, = K u ■ K u /k H . 

The up-type quark coupling-strength scale factor is mostly indirectly constrained through the gg —> H 
production channel, from the Higgs boson to top-quark coupling strength, with an additional weak direct 
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constraint from the qqlgg —> ttH production channel, while the down-type coupling strength is con¬ 
strained through the H —> bb, H —> tt and H —> pp decays as well as weakly through the bb —> H 
production mode and the 6-quark loop in the gg —> H production mode. 

The fit results for the parameters of interest in this benchmark model, when profiling the other parameters, 
are: 


A du e [-1.08, -0.81] U [0.75,1.04] (68% CL) 

A Vu = 0.92^ 

Kuu - l-25_ 033 . 

Near the SM prediction of A du = 1, the best-fit value is A du = 0.90 1 j i. All parameters are measured 
to be consistent with then - SM expectation and the three-dimensional compatibility of the SM hypothesis 
with the best-fit point is 51%. 

The likelihood curves corresponding to these measurements are shown in Figure 15. The likelihood 
curve of Figure 15a is nearly symmetric around A du = 0 as the model is almost insensitive to the relative 
sign of k u and K d . The interference of contributions from the 6-quark and /-quark loops in the gg —> H 
production induces an observed asymmetry of about 0.6cr (no significant asymmetry is expected with the 
present sensitivity). The profile likelihood ratio value at A du = 0 provides 4.5cr evidence of the coupling 
of the Fliggs boson to down-type fermions, mostly coming from the H —» rr measurement and to a 
lesser extent from the H —> bb measurement. Vanishing coupling strengths of the Fliggs boson to up-type 
fermions {k uu = 0) and vector bosons (Ay n - 0) are excluded at a level of > 5cr. 


5.3.2. Probing the quark and lepton symmetry 

Extensions of the SM can also contain different coupling strengths of the Fliggs boson to leptons and 
quarks, notably some variants of Two-Higgs-Doublet Models. In this benchmark model the ratio A/ q 
of coupling-strength scale factors to leptons and quarks is probed, while vector boson coupling-strength 
scale factors are assumed to be unified and equal to Ky. The indices l, q stand for all leptons and quarks, 
respectively. The free parameters are: 


A (q - K(/K q 
Avq = Ky/Kq 
Kqq - Kq • K q /K H . 

The lepton coupling strength is constrained through the H —> rr and H —> /rq decays. The fit results for 
the parameters of interest of this benchmark model, when profiling the other parameters, are: 

Af q e [-1.34, -0.94] U [0.94,1.34] (68% CL) 

Avq = 1.03!j>;H 
Kqq = 1-03 ^20- 

Near the SM prediction of A( q = 1, the best-fit value is Ac q = 1. 12 All parameters are measured to be 
consistent with their SM expectation and the three-dimensional compatibility of the SM hypothesis with 
the best-fit point is 53%. 
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(c) 


Figure 15: Results of fits for the benchmark model described in Section 5.3.1 that probes the ratio of scale factors 
between down- and up-type fermions: profile likelihood ratios as functions of the coupling-strength scale factor 
ratios (a) Ad„ ( Ay u and k uu are profiled), (b) Ay u ( Ad u and k uu are profiled), and (c) the overall scale factor k uu (Ad„ 
and Ay„ are profiled). The dashed curves show the SM expectations. The red (green) horizontal line indicates the 
value on the profile likelihood ratio corresponding to a 68% (95%) confidence interval for the parameter of interest, 
assuming the asymptotic^ 2 distribution for the test statistic. 


Figure 16 shows the likelihood curves corresponding to the fit results for this benchmark. Similar to the 
model of Section 5.3.1, the likelihood curve is nearly symmetric around Ae q = 0. A vanishing coupling 
strength of the Higgs boson to leptons, i.e. A( q = 0, is excluded at the ~ 4.4<x level due to the H —> tt 
measurement. The profile likelihood ratio values at K qq = 0 and Ay q = 0 provide strong confirmation of 
Higgs boson couplings to quarks and vector bosons with both significances of > 5<x. 


5.4. Probing beyond the SM contributions in loops and decays 

In this section, contributions from new particles either in loops or in new final states are probed. For 
the H—> yy, H —> Zy and gg —> H vertices, effective scale factors K y , Kz y and k (j are introduced that 
allow for extra contributions from new particles. These effective scale factors are defined to be positive as 
there is by construction no sensitivity to the sign of these coupling strengths. The potential new particles 
contributing to these vertex loops may or may not contribute to the total width of the observed state 
through direct invisible or undetected decays. In the latter case the total width is parameterised in terms 
of the additional branching ratio BRi., u . into invisible or undetected particles. 
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Figure 16: Results of fits for the benchmark model described in Section 5.3.2 that probes the symmetry between 
quarks and leptons: profile likelihood ratios as functions of the coupling-strength scale factor ratios (a) Af q (Ay q and 
K qq are profiled), (b) Ay q (Ac q and K qq are profiled), and (c) the overall scale factor K qq (A[ q and Ay q are profiled). The 
dashed curves show the SM expectations. The red (green) horizontal line indicates the value of the profile likelihood 
ratio corresponding to a 68% (95%) confidence interval for the parameter of interest, assuming the asymptotic x 2 
distribution for the test statistic. 


5.4.1. Probing BSM contributions in loop vertices only 

In the first benchmark model of this section, BSM contributions can modify the loop coupling strengths 
from their SM prediction, but it is assumed that there are no extra contributions to the total width caused 
by non-SM particles. Furthermore, all coupling-strength scale factors of known SM particles are assumed 
to be as predicted by the SM, i.e. kw = «z - K t - Kb - k t - k m - 1. The free parameters are thus K g , K y 
and Kz y . 

Figure 17a shows the results of fits for this benchmark scenario and the best-fit values and uncertainties, 
when profiling the other parameters. The effective coupling-strength scale factors k cj and K y are measured 
to be consistent with the SM expectation, whereas a limit is set on the effective coupling-strength scale 
factor Kzy Figure 17b shows the two-dimensional likelihood contour for n g vs. K y , where Kz y is profiled. 
The three-dimensional compatibility of the SM hypothesis with the best-fit point is 69%. 


5.4.2. Probing BSM contributions in loop vertices and to the total width 

The second benchmark model of this section removes the assumption of no invisible or undetected Higgs 
boson decays, introducing BRj u as additional model parameter. The free parameters of this benchmark 
model are thus K g , K y , Kz y and BR L U . The coupling-strength scale factors of known SM particles are 
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Figure 17: Results of fits for the benchmark model that probes for contributions from non-SM particles in the 
H --> yy, H —> Zy and gg —> // loops, assuming no extra contributions to the total width: (a) overview of fitted 
parameters, where the inner and outer bars correspond to 68% CL and 95% CL intervals, and (b) results of the 
two-dimensional fit to K y and k cj , including 68% and 95% CL contours (k/ 7 is profiled). 
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Figure 18: Results of fits for benchmark models that probe for contributions from non-SM particles in the H —> yy, 
H —> Zy and gg —> // loops, while allowing for potential extra contributions to the total width: (a) overview of 
fitted parameters. The inner and outer bars correspond to 68% CL and 95% CL intervals. The confidence intervals 
for BR; , u are estimated with respect to the physical bounds as described in the text, (b) Profile likelihood ratio 
as a function of the branching fraction BR; u to invisible or undetected decay modes (k 7 , k cj and Ky y are profiled). 
The red (green) horizontal line indicates the value of the profile likelihood ratio corresponding to a 68% (95%) 
confidence interval for the parameter of interest, assuming the asymptotic^ 2 distribution for the test statistic. 
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still assumed to be at their SM values of 1. Due to this assumption, the parameterisation of Higgs boson 
channels that do not involve a loop process, e.g. VBF production of H —> WW* and associated production 
of H —> bb , depends only on BRi u in this model, and not on Kg, K y or Kz y , and can hence constrain 
BRi u from the data. Thus no additional constraints, beyond those introduced in the benchmark model of 
Section 5.2.2, are necessary in this model. 

The results of fits to this benchmark model are shown in Fig. 18, along with the uncertainty on the total 
width that this model allows, obtained in the same fashion as for the previous benchmark models. The 
effective coupling-strength scale factors k cj and K y are measured to be consistent with the SM expectation, 
whereas limits are set on the effective coupling-strength scale factor Ky y and the branching fraction BR L L1 . 
By using the physical constraint BRi u > 0, the observed 95% CL upper limit is BR IU < 0.27 compared 
with the expected limit of BRj u < 0.37 under the SM hypothesis. The four-dimensional compatibility of 
the SM hypothesis with the best-fit point is 74%. The best-fit values of the model parameters of interest 
and their uncertainties, when profiling the other parameters, are 


*g = 1 • 12 -al! [^ao8( stat -) -ao5( s y st -) -o:o7( theo -)| 

K y = 1.00 + 0.12 [!aii(stat.) !ao 5 (syst-) !ao3(theo.)] 


In a variant of the fit where no limits are imposed on BR L ll its best-fit value is 


BRi .u. = -0-15 ±o; 22 [-ai7( sta t-)!o:ii(syst.)!o:o7(theo.)], 

corresponding to the likelihood curve shown in Fig. 18b. Without the condition BRi., u . > 0, the best- 
fit value of BR; -u assumes a small (unphysical) negative value that is consistent with zero within the 
uncertainty. 

As the choice of free parameters in this model gives extra degrees of freedom to ggF production and 
H—> yy and H —> Zy decays, the most precise measurements based on ggF production or H-* yy decays 
(see Fig. 2) do not give a sizeable contribution to the determination of BRj u . Instead BR IU is mostly 
constrained by channels sensitive to VBF and VH production, as the tree-level couplings involved in these 
production modes are fixed to their SM values within this model. The upward uncertainty on T///T™ is 
notably increased with respect to that of the model in Section 5.4.1 due to the removing the constraint on 
BRi., u ., whereas the downward uncertainty is identical due to the condition that BR, U > 0. 


5.4.3. Probing BSM contributions in loop vertices and to the total width allowing modified 
couplings to SM particles 

The last benchmark model of this section removes the assumption of SM couplings of the Higgs boson for 
non-loop vertices used so far in this section, re-introducing the coupling-strength scale factors Kp and Ky 
defined in Section 5.2.1 to allow deviations of the coupling strength of the Higgs boson to fermions and 
gauge bosons, respectively. As the expression for kh is no longer strongly constrained due to the newly 
introduced degrees of freedom, the upper limit on Yh is no longer bounded, and an additional constraint 
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on the total Higgs boson width must be introduced. Similar to the model of Section 5.2.2 the two choices 
of the constraints on the total width discussed in Section 5.1 are studied: Ky < 1 and K on = The free 
parameters of this model are Kp, Ky, k ( j , k 7 , kz 7 and BR; u . 

Figure 19 shows the best-fit values and their uncertainties. The coupling-strength scale factors k cj , k 7 , 
Ky and Kp are measured to be consistent with their SM expectation, while limits are set on the coupling- 
strength scale factor Kz y and the branching fraction BR L L1 to invisible or undetected decays. By using the 
physical constraint BRj u > 0, the 95% CL upper limit is BRj u < 0.27 (0.54) when applying the con¬ 
straint Ky < 1 (/c on = K 0 ff). The expected limit in case of the SM hypothesis is BRj u < 0.39 (0.72). The 
six-dimensional compatibility of the SM hypothesis with the best-fit point is 96% (64%) when applying 
the Ky < 1 (/Con = *"off) constraint, respectively. The uncertainty on Y ///T^ M is significantly increased com¬ 
pared with models in Sections 5.4.1 and 5.4.2 due to the further relaxed coupling constraints, in particular 
both the 68% and 95% CL intervals of T^/T™ extend below 1. 


5.5. Generic models 

In the benchmark models studied in Sections 5.2, 5.3 and 5.4, specific aspects of the Higgs sector are 
tested by combining coupling-strength scale factors into a minimum number of parameters under certain 
assumptions, thereby maximising the sensitivity to the scenarios under study. In generic models the scale 
factors for the coupling strengths to W, Z, t, b, r and n are treated independently, while for the loop 
vertices and the total width Th, either the SM particle content is assumed (Section 5.5.1) or no such 
assumption is made (Sections 5.5.2 and 5.5.3). 


5.5.1. Generic model 1: no new particles in loops and in decays 

In the first generic benchmark model all coupling-strength scale factors to SM particles, relevant to the 
measured modes, are fitted independently. The free parameters are: kw, kz, K t , Kb, k t , and K fJ . It is 
assumed that only SM particles contribute to Higgs boson vertices involving loops, and modifications 
of the coupling-strength scale factors for fermions and vector bosons are propagated through the loop 
calculations. No invisible or undetected Higgs boson decays are assumed to exist. Only the W coupling- 
strength scale factor is assumed to be positive without loss of generality: due to interference terms, the fit 
is sensitive to the relative sign of the W and t couplings (through the iH, H— »yy, H —» Zy processes) and 
the relative sign of the Z and t coupling (through the gg —> ZH process), providing indirect sensitivity 
to the relative sign of the W and Z coupling. Furthermore, the model has some sensitivity to the relative 
sign of the t and b coupling (through the ggF process). 

Figure 20 summarises the results of the fits for this benchmark scenario. All measured coupling-strength 
scale factors in this generic model are found to be compatible with their SM expectation, and the six¬ 
dimensional compatibility of the SM hypothesis with the best-fit point is 57%. Illustrative likelihoods 
of the measurements summarised in Fig. 20 arc shown in Fig. 21. As shown in Figs. 21a and 21b, the 
negative solution of K t is strongly disfavoured at 3.1cr (2.9cr expected), while the negative minimum of /q, 
is slightly disfavoured at 0.5cr (no sensitivity expected). 

For the measurements in this generic model, it should be noted that the low fitted value of k/, causes a 
reduction of the total width Th by about 30% compared to the SM expectation (see Table 9), which in 
turn induces a reduction of all other k- values by about 20%. 
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Figure 19: Results of fits for benchmark models that probe for contributions from non-SM particles in the II—> yy, 
gg —> H and H —> Zy loops, with free gauge and fermion coupling-strength scale factors k v , k f , while allowing 
for potential extra contributions to the total width. The estimated values of each parameter under the constraint 
Ky < 1 or K on = k 0 ff are shown with markers in the shape of a box or a circle, respectively. The inner and outer bars 
correspond to 68% CL and 95% CL intervals. The confidence intervals of BR; u and, in the benchmark model with 
the constraint «y < 1, also k v , are estimated with respect to their physical constraints as described in the text. The 
numerical values of the fit under the constraint Ky < 1 are shown on the left. Values for the alternative A: on = /c 0 ff 
constraint are also shown (in a reduced font size due to space constraints). 
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Figure 20: Overview of best-fit values of parameters with 68% and 95% CL intervals for generic model 1 (see text). 
In this model only SM particles are considered in loops and no invisible or undetected Higgs boson decays are 
allowed. The sign of k w is assumed to be positive, as indicated by the hatched area, without loss of generality. The 
inner and outer bars correspond to 68% CL and 95% CL intervals. 


Figure 22 shows the results of the fit for generic model 1 as reduced coupling-strength scale factors 


Uvj = 


gvj ,—iiivj 

k V j— = y/Kvj—^- 


( 10 ) 


for weak bosons with a mass my, where gyj is the absolute Higgs boson coupling strength, v is the vacuum 
expectation value of the Higgs field and 


gp,i mp,i 

ilFJ - Kfj—z = Kf i - 

V2 V 


( 11 ) 


for fermions as a function of the particle mass mp, assuming a SM Higgs boson with a mass of 125.36 GeV. 
For the ft-quark mass in Fig. 22 the MS running mass evaluated at a scale of 125.36 GeV is assumed. 
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(a) (b) 
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(d) 


Figure 21: Results of fits for generic model 1 (see text): profile likelihood ratios as functions of the coupling- 
strength scale factors (a) K t , (b) Kb, (c) k w , and (d) k z . For each measurement, the other coupling-strength scale 
factors are profiled. The kinks in the curves of (a) and (c) are caused by transitions in solutions chosen by the profile 
likelihood for the relative sign between profiled couplings. The dashed curves show the SM expectations. The red 
(green) horizontal line indicates the value of the profile likelihood ratio corresponding to a 68% (95%) confidence 
interval for the parameter of interest, assuming the asymptotic;^ 2 distribution for the test statistic. 


5.5.2. Generic model 2: allow new particles in loops and in decay 

In the second generic benchmark model the six free parameters from the first generic model are retained 
but the assumptions on the absence of BSM contributions in loops and to the total width are dropped. 
Effective coupling-strength scale factors for loop vertices are introduced, and optionally a branching ratio 
BRj u to new non-SM decays that might yield invisible or undetected final states is introduced, resulting 
in a total of 9 (10) free parameters. In the variant where BR, u is not fixed to zero, either the constraint 
Ky < 1 is imposed, or the constraint on the total width from off-shell measurements is included. 

Figure 23 summarises the results of the fits for this benchmark scenario. The numerical results are shown 
in Table 11. As an illustration of contributions from different sources, the uncertainty components are 
shown for the case of BRj u = 0. All fundamental coupling-strength scale factors, as well as the loop¬ 
coupling scale factors k cj and K y are measured to be compatible with their SM expectation under all 
explored assumptions, while limits are set on the loop-coupling scale factor Ky y and the fraction of Higgs 
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Figure 22: Fit results for the reduced coupling-strength scale factors yyy = y] K v,i ^7 = yfKyi'-f- for weak bosons 
and y Fi = k F j for fermions as a function of the particle mass, assuming a SM Higgs boson with a mass 

of 125.36 GeV. The dashed line indicates the predicted mass dependence for the SM Higgs boson. 


boson decays to invisible or undetected decays. When imposing the physical constraint BR; u > 0 in 
the inference on BRj u , the 95% CL upper limit is BRj u < 0.49 (0.68) under the constraint Ky < I 
(/ton = /toff) on the Fliggs boson total width. The nine-dimensional compatibility of the SM hypothesis 
with the best-fit point is 73% when BRj u is fixed to zero. The compatibilities for the tits with the 
conditions Ky < I and /c on = /c 0 ir imposed are 80% and 57%, respectively. 

Similar to the results of the benchmark model in Section 5.2.2 the upper bound of the 68% CL interval for 
the scenario /c on = should be considered to be only approximate due to deviations of the test-statistic 
distribution from its asymptotic form. The deviation of the asymptotic distribution was shown to be 
negligible for off-shell signal strengths corresponding to the upper end of the 95% asymptotic confidence 
interval (Table 11). 

Also shown in Fig 23 are the resulting ranges of the total width of the Fliggs boson, expressed as the ratio 
r H /T s ™. These estimates are obtained from alternative parameterisations of these benchmark models, 
where the effective coupling-strength scale factor K g is replaced by the expression that results from solving 
Eq. (8) for K g , introducing r///r^ M as a parameter of the model. The figure shows that the upper bound 
on the Fliggs boson width from the assumption K cin = k 0 s is substantially weaker than the bound from the 
assumption Ky < 1. These results on T#/T^ M represent the most model-independent measurements of the 
Fliggs boson total width presented in this paper. 

Figure 24 shows profile likelihood ratios as a function of selected coupling-strength scale factors. In 
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Figure 23: Results of fits for generic model 2 (see text): the estimated values of each parameter under the constraint 
Ky < 1, A'on = Kcff or BRj u = 0 are shown with markers in the shape of a box, a circle, or a diamond, respectively. 
The hatched area indicates regions that are outside the defined parameter boundaries. The inner and outer bars 
correspond to 68% CL and 95% CL intervals. The confidence intervals of BR; u and, in the benchmark model with 
the constraints kw < 1 and \k-/\ < 1, also k w and Ky, are estimated with respect to their physical bounds as described 
in the text. Numerical results are shown in Table 11. 
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(C) (d) 

Figure 24: Results of fits for generic model 2 (see text): profile likelihood ratios as functions of the coupling- 
strength scale factors (a) k,, (b) Kb, (c) k w , and (d) k z . For each measurement, the other coupling-strength scale 
factors are profiled. The red (green) horizontal line indicates the value of the profile likelihood ratio corresponding 
to a 68% (95%) confidence interval for the parameter of interest, assuming the asymptotic x 2 distribution for the 
test statistic. 


Fig. 24a, the negative minimum of K t is shown to be disfavoured at l.Ocr. The minimum corresponding to 
the positive solution is found at K t = 1.28 The sensitivity to disfavour the negative solution of k, is 
reduced with respect to generic model 1 as the interference in loop couplings can no longer be exploited 
because effective coupling-strength scale factors were introduced. The observed residual sensitivity to the 
sign of Kf is exclusively due to the tree-level interference effect of the tH background in the ttH channel. 

The power of individual loop processes to measure the magnitude of K t and resolve the sign of k, relative 
to Kw is illustrated in more detail in Fig. 25. The blue curve shows the profile likelihood ratio as a 
function of K t for a model with the least sensitivity to the sign of k,: all loop processes are described with 
effective coupling parameters, including the gg —» ZH loop process. Subsequently the red, green and 
orange curves represent the profile likelihood ratios for models that incrementally include information 
from loop processes by resolving the gg —» ZH, ggF and H —> yy, Zy loop processes into their expected 
SM content. Here the red curve corresponds to the configuration of generic model 2, and the orange 
curve corresponds to the configuration of generic model 1. As expected, resolving gg —> ZH process 
adds little information on K t . Additionally resolving the ggF loop process into its SM content greatly 
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improves the precision on k, (green curve), but reduces the sensitivity to the relative sign of k, and kw- 
This reduction happens because on one hand the ggF process yields no new information on this relative 
sign, as it is dominated by t-b interference, and on the other hand because it decreases the observed 
magnitude of k, to a more SM-compatible level, thereby reducing the sensitivity of the iH process to the 
relative sign. Further resolving the H—> yy and H —> Zy loop processes, which are dominated by W-t 
interference, greatly improves the measurement of the relative sign of kw and k, (orange curve), but does 
not significantly contribute to the precision of the magnitude of k,. 



Figure 25: Profile likelihood ratio as a function of k, for models with and without resolved loop processes: shown 
are measurements of k, with no loop processes resolved (blue), only gg —» ZH resolved (red, generic model 2), 
gg —> H additionally resolved (green), and H-* yy and H —> Zy additionally resolved (orange, generic model 
1). The dashed blue and orange curves correspond to the expected sensitivity for the no-loop and all-loop models. 
All profile likelihood curves are drawn for the full range of K t , however some curves are partially obscured when 
overlapping with another nearly identical curve. The red (green) horizontal line indicates the value of the profile 
likelihood ratio corresponding to a 68% (95%) confidence interval for the parameter of interest, assuming the 
asymptotic^ 2 distribution for the test statistic. 


5.5.3. Generic model 3: allow new particles in loops, no assumptions on the total width 

In the final benchmark model of this section, the six absolute coupling-strength scale factors and three 
effective loop-coupling scale factors of generic model 2 are expressed as ratios of scale factors that can 
be measured independent of any assumptions on the Higgs boson total width. The free parameters are 
chosen as: 
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K gZ 

A Zg 

= K g ■ K Z /k H 

= KzlKg 

4rZ 

= k t /k z 

A wz 

= K W /Kz 

T//Z 

= Kfj/Kz 

4-tg 

= Kfl K CJ 

AyZ 

= Ky/Kz 

A bZ 

= K b /Kz 

T(Zy)Z 

= KZy / Kz 


Figure 26 shows the full set of results obtained from the fit to this benchmark model. The fitted values and 
their uncertainties are also shown in Table 12. As the loop-induced processes arc expressed by effective 
coupling-strength scale factors, there is little sensitivity to the relative sign of coupling-strength scale 
factors due to tH and gg —> ZH processes only. Flence only positive values for all /r- factors except k, 
are shown without loss of generality. The parameter K g z, Azg, Awz, d tg, d tz , d r z and A y z are all measured 
to be compatible with their SM expectation, while limits are set on the parameters A g z and A(z y )z■ The 
nine-dimensional compatibility of the SM hypothesis with the best-fit point is 73%. 

Table 12: Numerical results of the fits for generic model 3: measurements of ratios of coupling-strength scale factors 
in which assumptions on the Higgs boson total width cancel. These results are also shown in Fig. 26. Shown in 
square brackets are uncertainty components from different sources. For Awz and A tg , the uncertainty breakdowns 
are provided for the preferred positive solutions. 


Parameter Measurement Uncertainty breakdown 


KgZ 

= 

1.18 + 0.16 

!al|( s tat.) 

-no4( s y st -) 

:^(theo.) 

A Zc, 

= 

1 09 +0 - 26 

1 - u -0.22 

-o'2C)( Stat ') 

t| 2 (syst.) 

to6( theo -). 

A WZ 

6 

[-1.04,-0.81] U [0.80,1.06] 

— o’l i (stat.) 

^“(syst.) 

-(I02( theO -) 

A tg 

€ 

[-1.70,-1.07] U [1.03,1.73] 

:« ; g( s tat.) 

loi^syst.) 

-o'of|( the0 ') 

AhZ 

= 

0.60 ± 0.27 

-o'i9( stat ') 

^(syst.) 

-o'o3 ( the °') 

TrZ 

= 

0 qo+0.23 
u ' -0.19 

-0+6( Stat -) 

-o:o9( s y st -) 

-o'o^theo-) 

ITpzl 

< 

2.3 (95% CL) 




TyZ 

= 

0.90 + 0.15 

^(sfat.) 

Iao4( s y st -) 

-o'o3 ( theo ') 

|d(Zy)zl 

< 

3.2 (95% CL) 





The parameter Awz - K w! K z in this model is of particular interest: identical coupling-strength scale 
factors for the W and Z bosons are required within tight bounds by the SU(2) custodial symmetry and the 
p parameter measurements at LEP and at the Tevatron [120]. This custodial constraint is directly probed in 
the Fliggs sector through the parameter Awz- The measured ratio Awz is in part directly constrained by the 
decays in the H—> WW*—>Mv and H—> ZZ*—> 4f channels and the WH and ZH production processes. It 
is also indirectly constrained by the VBF production process, which in the SM is 74% W fusion-mediated 
and 26% Z fusion-mediated (see Table 9). Figure 27a shows the profile likelihood ratio as a function of the 
coupling-strength scale factor ratio Awz- Due to the interference terms, the fit is sensitive to the relative 
sign of the W and t coupling (tH) and the relative sign of the Z and t coupling (gg —> ZH ), providing 
indirect sensitivity to the sign of Awz- The negative solution is disfavoured at 0.5cr (0.3cr expected). The 
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minimum corresponding to the positive solution is found at Awz = 0.92^[j‘ ] ^, in excellent agreement with 
the prediction of SU(2) custodial symmetry. 

Also shown in Figs. 27b, 27c are the ratios A y z and A lg . The ratio A y z is sensitive to new charged particles 
contributing to the H—»yy loop in comparison to H —* ZZ* decays. Similarly, the ratio A, g is sensitive 
to new coloured particles contributing through the gg —» H loop as compared to ttH. The minimum 
corresponding to the positive solution is found at A tg - 1.38 ± 0.35. Both arc observed to be compatible 
with the SM expectation. 

The fit in the third generic benchmark model uses only the basic assumptions, as stated at the beginning of 
this section, and hence represents the most model-independent determination of coupling-strength scale 
factors that is currently possible. 
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K gZ = 1.18 ±0.16 

^= 1 - 09^11 

A wz e[- 1.04,-0.81] 
u[0.80,1.06] 
A tg e [-1.70,-1.07] 
u[1.03,1.73] 

A bZ = 0.60 ±0.27 

A rZ = 0.99 !q;?9 
(95 %CL) A^z <2.3 
A y z = 0.90 ± 0.1 5 
(95 %CL) A {Zy)Z < 3.2 
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Parameter value 


Figure 26: Results of fits for generic model 3 (see text): allowing deviations in vertex loop-coupling scale factors and 
in the total width. Overview of best-fit values of parameters, where the inner and outer bars correspond to 68% CL 
and 95% CL intervals. The hatched areas indicate regions that are outside the defined parameter boundaries. 
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Figure 27: Results of fits for generic model 3 (see text): profile likelihood ratios as functions of the coupling- 
strength scale factor ratios (a) Awz, (b) A tg and (c) A y z- In all cases, the other parameters are profiled. The dashed 
curves show the SM expectations. The red (green) horizontal line indicates the cutoff value of the profile likelihood 
ratio corresponding to a 68% (95%) confidence interval for the parameter of interest, assuming the asymptotic % 2 
distribution for the test statistic. 
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6. Conclusion 


The Higgs boson production and decay properties are studied using proton-proton collision data collected 
by the ATLAS experiment at the Large Hadron Collider corresponding to integrated luminosities of up 
to 4.7 fb -1 at sjs = 7 TeV and 20.3 fb _I at = 8 TeV. The study combines specific analyses of the 
H — > yy, ZZ*, WW*, Zy, bb, tt and p// decay channels, as well as searches for ttH production and 
measurements of off-shell Higgs boson production. It significantly extends a previous combination of 
the H —» yy, ZZ* and WW* decays [23]. In particular, the addition of the fermionic decays of the Higgs 
boson in the combinations allows for direct tests of the Yukawa interactions of the Higgs boson with 
fermions. 

The measured Higgs boson signal yields are compared with the SM expectations at the fixed Higgs bo¬ 
son mass of mu = 125.36 GeV. The combined yield relative to its SM prediction is determined to be 
1.18 ± 0.10 (stat.) ± 0.07 (syst.) (theo.). The combined analysis provides unequivocal confirmation 
of gluon fusion production of the Higgs boson with a significance exceeding 5cr and strong evidence of 
vector-boson fusion production with a significance of 4.3cr. Furthermore, it supports the SM predictions 
of Higgs boson production in association with a vector boson or a pair of top quarks. Values for the total 
cross sections can be obtained from the signal strength of each production process within the uncertain¬ 
ties related to the modelling of Higgs boson production and decay kinematics and assuming SM decay 
branching ratios. The total cross sections at Vi = 7 and 8 TeV are 22.1 ^ ^ (stat.) ^'3 (syst.) (theo.) pb 
and 27.7 ± 3.0 (stat.) (syst.) (theo.) pb, respectively. 

The observed Higgs boson production and decay rates are also interpreted in a leading-order coupling 
framework, exploring a wide range of benchmark coupling models both with and without assumptions 
about the Higgs boson width and the SM particle content of loop processes. Higgs boson couplings to up- 
type fermions and vector bosons are found with both significances above 5 cr and to down-type fermions 
with a significance of 4.5cr, under the assumption of unified coupling scale factors, one for each type of 
particles. In a different model with separate unified coupling scale factors for leptons, quarks and vector 
bosons, Higgs boson couplings to leptons are found with a significance of ~ 4.4cr. 

The Higgs boson coupling strengths to fermions and bosons are measured with a precision of ±16% and 
±7% respectively, when assuming the SM Higgs boson width, and are observed to be compatible with 
the SM expectations. Coupling strengths of loop processes are measured with a precision of ±12% when 
assuming the SM expectations for non-loop Higgs boson coupling strengths and the Higgs boson total 
width, increasing to about ±20% when these assumptions are removed. No significant deviations from 
the SM expectations of Higgs boson coupling strengths in loop processes are observed. 

Measurements of coupling strengths to /j, r leptons, b, t quarks and W, Z bosons, or ratios of these coup¬ 
ling strengths are presented in the context of generic Higgs boson coupling models. They can constrain 
the ratio of W and Z coupling strengths, a probe of custodial symmetry, with a precision of ±13%. For 
benchmark models that measure absolute coupling strengths, a variety of physics-motivated constraints on 
the Higgs boson total width have been explored. The measured Higgs boson coupling strengths and their 
precision are found to depend only weakly on the choice of these constraints. A third generic benchmark 
model uses only the most basic assumptions and hence represents the most model-independent determin¬ 
ation of the coupling strength scale factors that is currently possible. In this model ratios of couplings are 
constrained with a precision of 15 - 40%. 
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The /^-values expressing compatibility of the SM hypothesis with the best-fit point range between 29% 
and 99% for all considered benchmark models. The observed data arc thus very compatible with the SM 
expectation under a wide range of assumptions. 
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Appendix 

A. Alternative parameterisation of ratios of cross sections and of 
branching ratios 


An alternative to the parameterisation of Section 4.4 is to normalise the ratios of cross sections and of 
branching ratios to their SM values. Compared with Eq. (6), the yield of the production and decay 
i —* H —* f can be parameterised as 

cr• BR f = p{ X [a ■ BR/] sm - (/r^p • R i/esF • p f /ww) x [cr,- • BR/] SM . (12) 


Here R and p arc ratios of cross sections and branching ratios relative to their SM expectations, respect¬ 
ively: 


oy/crggp 

K i/ ggF - 7 - i — ana Pf/ww 

[o-i/o- ggF J SM 


BRy/BRww* 

[bR//BRww*] sm 


(13) 


The data are fitted with p w ^ , four ratios of production cross sections and one ratio of branching ratios 
for each decay channel other than the H —> WW* decay. The results shown in Table 13 are nearly 
identical to the best-fit values relative to their SM predictions shown in Table 7. The small differences 
are expected from the inclusion of additional nuisance parameters of the SM predictions and from the 
precision of the fits. One clear advantage of the parameterisation of Section 4.4 is that the results arc 
independent of the SM predictions and are, therefore, not affected by the theoretical uncertainties of the 
predictions. Consequently, the fitted values of the ratios of cross sections and of partial decay widths 
shown in Table 7 have significantly smaller theoretical uncertainties than their counterparts (R,/ gg F and 
Pr/ww *) in Table 13. The former is only affected by the theoretical uncertainties in the modelling of Higgs 
boson production whereas the latter suffer from both the modelling uncertainties and the uncertainties of 
the SM predictions. 
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Table 13: Best-fit values of gg —» H —» WW* signal strength p w ^\ ratios of cross sections /?,■/ gg p and of branching 
ratios pf/ww- All /?,• / gg F and Pf/ww are measured relative to their SM values for m H = 125.36 GeV from the 
combined analysis of the yfs — 7 and 8 TeV data. Shown in square brackets are uncertainty components: statistical 
(first), systematic (second) and signal theoretical (third) uncertainties. 


Parameter Best-fit value 


..WW 

^ggF 

1 1 5 

+0.28 

+0.18 

+0.12 

+0.171 


-0.24 

-0.18 

-0.11 

-0.12 J 

^VBF/ggF 

n qq 

+0.46 

+0.37 

+0.20 

+0.181 

\J.yy 

-0.33 

-0.29 

-0.12 

-0.101 

R WH/ ggF 


+ 1.06 

+0.87 

+0.49 

+0.34] 

1.4/ 

-0.74 

-0.65 

-0.32 

-0.15 | 

R ZH/ ggF 

0.60 

+ 1.39 
-0.66 

+0.99 

-0.60 

+0.93 

-0.25 

+0.30] 

-0.07] 

R ttH/ ggF 

1.81 

+1.10 

-0.81 

+0.79 

-0.64 

+0.61 

-0.48 

+0.46] 

-0.17] 

Pyy/WW 

0.97 

+0.32 

-0.25 

+0.26 

-0.22 

+0.15 

-0.10 

+0.101 
-0.06] 

PZZ’/WW 

1 0/1 

+0.42 

+0.37 

+0.18 

+0.07] 

1 .Z4 

-0.31 

-0.29 

-0.10 

-0.04] 

Ptt/WW* 

1.20 

+0.52 

-0.38 

+0.40 

-0.32 

+0.29 

-0.18 

+0.171 
-0.09] 

Pbb/WW* 

0.59 

+0.63 

-0.37 

+0.45 

-0.27 

+0.43 

-0.24 

+0.12] 
-0.05 1 
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